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ABSTRACT 


A lean combustion model for monolithic catalytic combustors is 
given. The model, consisting of several semi-global chemical reaction 
steps in the gas phase and on the surface, is capable of analyzing CO 
and THC emissions and combustor efficiency in both steady and unsteady 
states . 

In the steady-state model computation presented, the influence 
of operating and design parameters on the minimum combustor length 
is studied. Special attention is given to the effect of after-bed gas 
phase reaction space. Comparison with experimental data by Anderson 
indicates good agreement in the range of parameters covered. 

In the transient analysis, a quasisteady gas phase and a thermally- 
thin substrate are assumed. The combustor response delay is due to the 
substrate thermal inertia. Fast response is found to be favored by thin 
substrate, short catalytic bed length, high combustor inlet and final 
temperatures and, in most cases, small gas channel diameters. The 
calculated gas and substrate temperature time history at different axial 
positions provide an understanding of how the catalytic combustor responds 
to an upstream condition change. The computed results also suggest that 
the gas residence times in the catalytic bed and in the after-bed space 
are correlatable with the nondimensional combustor response time. When 
fast transient responses are required, both steady and unsteady studies 
have to be made to achieve a meaningful compromising design. 
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INTRODUCTION 


The operation and the design of catalytic combustors are limited on 
one end by the high temperature that the catalysts can tolerate and on 
the other end (low temperature) by the acceptable emission levels of CO 
and unburnt hydrocarbons . The aim of the combustion engineer is to design 
a compact and durable combustor capable of efficient and clean combustion. 
Catalytic combustor modeling can contribute to this objective by providing 
more understanding to the detail processes occurring in the reactor and to 
serve as a design guidance for minimizing the numbers of experimental tests 
needed . 

There are several steady-state catalytic combustor models in existence 
(Ablow and Wise, 1979, Cerkanowicz, et al, 1977, Kelley et al, 1977, Kendall 
et al, 1979) and they show a high degree of success in describing the events 
occurring in the combustor. In particular, the model by Kelley et at (1977) 
containing detailed gas phase chemistry for methane/air reactions can 
predict pollutant emissions. Detailed kinetics, unfortunately, are not 
yet available for most other practical fuels. Since emission characteris- 
tics are important considerations for catalytic combustor design, in the 
absence of more detailed information, models consisting of several key 
semi-global chemical steps seem appropriate. 

In addition to the emission characteristics, the transient behavior 
of catalytic combustors is expected to be different from that of a conven- 
tional gas turbine burner, due to the large thermal Inertia of the substrate. 
Non-steady operations are obviously Important to transportation engines 
since their power levels have to be changed frequently. Transients may 
also be a concern with stationary gas turbine applications, since the 



ignition/shutdown operation can produce excessive thermal stress or thermal 
shock if the combustor is not designed or operated properly and the substrate 
can fail as a result (DeCorso and Carl, 1979). 

Typical questions concerning transient operations include the combustor 
response time, the type of response, the unsteady substrate temperature 
history and the resultant thermal stress distribution. At the present time, 
both experiment and theory on transient catalytic combustion are lacking. 

A model is presented here attempting to analyze the problem theoretically.. 

In Chapter I of this report, the transient combustion model will be formulated 
and numerical solution will be presented as a parametric study. Since the 
transient computer code can also yield steady-state results, in Chapter II, 
a detailed computation of steady states is given and a comparison with 
Anderson's data is made. As far as the combustion models are concerned, the 
one in Chapter II has a three-step global gas-phase chemical reaction vs. 
a two-step one in Chapter II. The reason for using the three-step scheme 
is discussed in the beginning of Chapter II. 
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CHAPTER I TRANSIENT CATALYTIC COMBUSTOR MODEL 

I.l Nonsteady Combustion Model 

Consideration of Time Scales 

For a monolithic catalytic combustor with uniform cell distribution In 
the plane normal to the approaching flow direction. If the combustor wall Is 
well Insulated, then the study of the whole reactor can be reduced to the con- 
sideration of a single cell unit. Referring to Fig. I.l, each cell unit 
consists of an open gas channel, the associated substrate volume and an 
after -bed gas reaction space. 
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Figure I.l Schematic Drawing of One Cell Unit 
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Considering first the steady-state combustion, the relevant time 
scales are the gas residence time in the reactor, the heat and mass transfer 
time between the gas and the solid surface inside the channel and the gas- 
phase and surface reaction times. When a transient is caused by an upstream 
parameter variation, then the new time scales Involved Include the input 
characteristic time, the time for the temperature wave to reach the central 
plane of the substrate and the time for heating-up (or cooling-down) the 
substrate. 

Table LI defines these characteristic times and their estimated magni- 
tudes. It should be noted that some of the time scales vary along the 
channel of the reactor bed. For example, the heat/mass transfer time and 
the substrate heat-up time are shorter at the flow entrance region than 
those for far downstream, and the ones used in Table I are those based on 
fully-developed flow profiles. However, the estimation from Table I does 
serve the purpose of order-of-magnitude comparison. From this table we 
see that the longest time is the solid heat-up time, being of the order 
of seconds. If the time scale we are interested in is much greater than 
the gas residence time (~15 msec), then the gas phase processes (including 
heat/mass transfer and chemical reactions) can be regarded as in a quasi- 
steady state. If furthermore the substrate half-thickness is smaller 
than about .2 mm, then the substrate temperature distribution in the 
direction normal to the channel axis can be regarded as uniform at any 
given time. In this limiting case then, a "quasi-steady-gas-phase and thermally' 
thin-substrate" model is applicable with the only transient process being 
the substrate heating-up or cooling-down. In the following section, this 
model will be described in mathematical form. 

Quaslstaady Gas Phase 

Quasisteady gas phase implies that the differential equations for the 
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Table I.l Estimate of Transient Time Scales in Catalytic Monolith Combustor 


Time 

Definition 

Estimated Magnitude 

Gas residence time 
reactor bed 

1 

in 1 

1 

(L*/u*^f) (A*/A*) 

1-15 msec 

Gas residence time 
after-bed space 

in ! 

! 

1 

L* /u* ^ 
as ref 

i 

0-15 msec 

Heat/mass transfer 

time 

11 d*^ 

Nu^ Tr 4a*(0,0) 

.5 - 7.5 msec for 

between gas and solid 
surface in channel j 

d* < 3.6 mm 

Chemical reaction times 
(gas-phase and heteroge- 


Same order of magnitude 
as gas residence times 

neous) 




Time for temperature wave 
to reach the substrate 
center plane 

2 

d* 

1 s 

a* 

1 s 

1 

0.5 - 25 msec for 

d* < .2 mm 
s 

Substrate heat-up 
time 


1 , - A*p*C* 

, 1 1 s s s 

1 Nu^ TT k*(0,0) 

I 

0.5 sec - 20 sec 

Transient input 
time 


{ 

j 

variable, specified 

i 

\ 


I 
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gas-phase processes are the same as those for the steady state. Only the 
boundary conditions on the solid surface are different. There are several 
steady-state catalytic combustor models in existence [Cerkanowicz et al, 1977, 
Kelly et al, 1977, 1979] and they show a high degree of success in describing 
the events occurring in the combustor. In particular, the model by Kelley et al 
(1977) containing detailed chemistry for methane/air reactions can predict 
pollutant emissions. Detailed kinetics, unfortunately, are not yet available 
for most other practical fuels. Since emission characteristics are important 
considerations for catalytic combustor design, in the absence of more 
detailed Information, a model consisting of several key semi-global chemical 
steps is proposed. 

For fuel-lean catalytic combustion using nitrogen-free hydrocarbon 
fuels, the NO emission is negligible because of low flame temperatures. 

The only pollutants needed to be considered are unburnt hydrocarbons (UHC) 
and carbon monoxide (CO) . For this reason, the two semi-global chemical 
reactions used in the gas phase are: 


(A) C H + 1/2 (n + m/2) 0„ n CO + m/2 H„0 
n m / / 


d[C H ] 
n m 


J ^ - c* p*“l [C H ]“2[o ]“3 exp (-E*/RT*) 

d t* 1 ■■ n m 2 1 

( 1 ) 


(B) 


CO + 1/2 O 2 CO^ 


= _C* p*^l [CO]^2j-o^]e3[H^O]®5 x*^4 (_e*/RT*) 

( 2 ) 


In describing the flow in the reactor channel, we follow previous 
investigators in using the plug flow approximation. For transient cases, in 
general, the fluid properties are functions of x, the distance along the 
channel axis and t, the time. For example, density of the gas will be 
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expressed as p(x,t). The nondimensional forms of these equations are 
summarized in the following (supscript * denotes dimensional quantity 
and without * denotes nondimensional quantities) . 


Continuity : 

pu = p(o,t) u(o,t) 

(3) 

Momentum : 

p = p(o,t) 

(4) 

Energy: 

pu |1 + J^(x . x^) = qiYg^(0,0) »1 



+ q2\^(0,0) W2 (5) 


Species : 

Hydrocarbon: 


Pu + J ( ‘Jhc - 3^HC,S) = - Wj (6) 

3x 


Carbon monoxide: 
pu 


9 ^£ 9 . + j ( ^co - ^co,s) = ’"iC - ^^2 

-- U 2 n 


3x 


(7) 


Oxygen : 


= 1 - <K1- “^HC) ^CO ( 




( 8 ) 

(9) 


Equation of state: p = pT 

The reaction rates are given by 

V. . B, ^-(Ei/T) 


W2 = B2 p ^ p 


( 10 ) 




( 11 ) 


and V ,^^(0.0) j 3 ,hc<*.'>] 

'^HC (12) 

The definition of and B 2 can be found in the nomenclature. 
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In deriving Eqns (5-8), heat conduction and mass diffusion in the 
axial direction are neglected because the Peclet number based on typical 
gas velocity (>10 M/sec) is much greater than unity [Ablow and Wise, 1979]. 
The dimenpionless lateral heat and mass transfer coefficients J and Jp.., 
in Eqns (5-7) are derived from Nusselt number calculation in entrance 
flow in a tube of constant surface temperature [Kays, 1966] . Following 
Kelly et al (1977), a modification is made at the entrance point based on 
local stagnation flow estimation. The use of entrance flow transport 
properties, rather than fully-developed constant values, make a great 
difference in the temperature and species distributions in the flow entrance 
region. 

The above system of equations can also be applied to the after-bed 
space where only gas-phase reactions occur. This is done by putting 
and equal to zero in Eqns (5-7) and by properly changing the mass flux 
per unit area in Eq (3) due to area change. 


Unsteady Solid with Thin-Wall Substrate 

When the substrate wall is thermally thin, the solid temperature can 
be regarded as a function of axial distance and time only, i.e., T^(x,t). 
The energy balance of the solid, including the solid surface, results in 
the following equation: 




8T' 


A”p:C„ - h*S*(T*- Tp = S -[q^jg + q^J^] 




* * 


(13) 


S S S 

In Eq (13) the heat conduction in the axial direction in the solid is 
neglected in comparison with the heat transfer rate between the gas and 
solid across the gas channel. 
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Two semi-global catalytic surface reactions are assumed. They are 


(C) + <n + f) O 2 * nCOj + f HjO 


(14) 


with the corresponding surface reaction rate J* given by 

J* - C* tc B 1 
= ^ n m s 


and (D) CO + -j °2 


(15) 




with 


* _* r-^-i — Eix/RT 

Ju = Cl [CO] e ^ 


In surface reaction (C) , the hydrocarbon is oxidized to foirm'CO^ 
and H 2 O, not CO such as in gas-phase reaction (A). This is based on 
the present available experimental evidence in oxygen-rich systems, 
only CO^, not CO, is found on or close to the catalytic surface [Schwartz 
et al, 1971, Anderson, 1976 and Marteney, 1979]. 

Eq (13) is non dimensionalized by first defining a proper reference 
time scale t*. 


r* = 


Nu 


1 

TT 


* * * 
A p C 
.S „ .S 


k*(0,0) 


(16) 


and 


*, * 

t /t 


_ Nu , * 
r , = X k 

nk 


Nu k*(0,0) 


(17) 


Eq (13) becomes 
9T 


s 


9 t 


nk 


(T - I^) - BjY^^CO.O) qjP eKp(-E3/T^) 


( 18 ) 
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Neglecting transient accumulation on the surface, the surface 
hydrocarbon and carbon monoxide are given by 


‘’nk ^ ^HC " "^HC,s^ “ O B 3 P VhC.s ®^P(-E 3 /Tg) 

^k< yco- ^CO,s>- Tl^3 ^C0,s 


Initial and Boundary Conditions 

The system of equations (3-12, 18-20) neeis initial and boundary 
conditions. Initial conditions (at t = 0) ate required for T^, ( 

( ^ function of x and upstream boundary conditions (at x = 0) 

should be specified for ^co’ ^02’ ^ u as a function of 

time. These conditions depend on the types of transient input which can 
vary with engine designs and mode of operations. In the present paper, 
a simple start-up transient is investigated, i.e. catalytic combustor 
response to a step-wise fuel Injection. 

The substrate temperature at time t=0 is specified. This temperature 
can be the result of external heating of the substrate, say, by a torch. 
The fuel and air are then supplied at upstream. The model is intended to 
describe the events from t*=0 to the final steady state. Mathematically, 
the conditions used are 
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= specified 


Tg(x,0) 

^HC, = - “ 

T(0,t) = 1 

- 1 
= ° 

'^02(0»t) = 1 

p(0,t) = 1 

u(0,t) = 1 


( 21 ) 


Numerical Solution Procedure 

Although most dependent variables in this problem are function of 
both X and t which, in general, results in partial differential equations, 
the assumption of quasisteady gas phase and thennally thin substrate 
greatly simplifies the mathematical property of the system and results 
in a simpler solution procedure. 

The gas-phase differential equations (5-7) have only first derivative 
in X with t as a parameter, the solid equation (18) has only derivative in 
t with X as a parameter; therefore, they can be integrated at each time (t) 
or position (x) as ordinary differential equations using Runge-Kutta method. 
The procedure is started by first integrating Eqns (5-7) forward in x from 
upstream to the end of combustor for T, ^ 


using the solid surface 
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quantities specified by the initial conditions, then Eq (18) is integrated 

for one time step (At) to find T (x. At) and Eqns (19,20) are solved 

6 

algebraically to find ^t) and '^qq g(x. At). These surface 

quantities are substituted back to Eqns (5-7) to start another cycle of 
Integration. Steady-state and transient solutions are obtained using the 
same code. The computational time ranges from 10-50 seconds in a VAX-11 
computer. The numerical scheme is an efficient one and it is suitable for 
parametric studies. 

1.2 Computed Results 
Steady Stat e 

Selected steady-state results will be presented first to facilitate 
the discussion on the transient response. Specifically, the relative 
importance of the catalytic bed length and the af ter-bed-gas-phase-reaction 
distance will be shown. In previous theoretical models, emphasis was placed 
on the combustion inside the catalytic bed only. The contribution of gas- 
phase reactions in the downstream space after the bed has been demonstrated 
in the experiments by T'ien and Anderson (1979) and Anderson (1980). 
Utilization of the pure gas-phase reactions in the downstream of the monolith 
to replace a portion of the catalytic bed will be shown to be crucial in 
reducing the catalytic combustor transient response time. 

Table 1. 2 lists the values of chemical kinetic constants for reactions 
A to D used in this ahap ter . The hydrocarbon fuel chosen is propane. Semi- 
global kinetics for hydrocarbon and CO oxidation are from Edelman (1969) and 
Dryer and Classman (1973) respectively, with the pre-exponential factors 
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Table 1.2 Chemical Kinetic Constants for Chapter I 


Reaction (A) 


Reaction (B) 


C* = 1,5 X 10^ C* = 0.71 X 10^^ 


= 0.3 


6 ^ = 0 


= 0.5 


^2 = ' 


“3 = ^ 


33 = 0.25 




^ = 0 


E* = 24 
(kcal/mole) 


= 0.5 


= 40 

(kcal/mole) 


Reaction (C) 

C* = 2.5x 10^ 

E* = 10 
(kcal/mole) 


Reaction (D) 

C* = 10^ 

4 

E* = 17.8 

4 

(kcal/mole) 
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adjusted. The activation energies for reactions C and D are obtained from 
the data of Marteney (1979) and Kuo and Morgan (1971) respectively with 
adjusted pre-exponential factors. It should be noted that the surface reaction 
rates contain information on catalyst type and loading density as well as 
the washcoat properties. 

The design parameters considered are d*, channel hydraulic diameter, 
A*/A^, the percentage of open area (or equivalently, A*/A*, the ratio of 
open-to-close areas), L*, the catalytic bed length and L* , after-bed gas-phase 
reaction length. The minimum catalytic combustor length required to reach 
the emission goal in steady-state operations is denoted using the subscript 
EG. The emission goal is chosen to be 1.64 g/kg of fuel for UHC and 13.6 
g/kg of fuel for CO [Anderson, 1977]. The operating parameters include p*, 
the pressure, combustor inlet temperature, T*^, adiabatic flame 

temperature (or equivalently the fuel/air equivalence ratio (()) and the 

reference velocity of the approach gas upstream of the reactor bed. 

Fig. 1.2 gives the profiles for one sample calculation. In this case, 
the propane fuel is oxidized quickly on the catalyst surface upon entering 
into the reactor channel and resulting in a high surface temperature, the 
gas temperature is raised by heat transfer from the substrate, gas-phase 
oxidation of propane is accelerated to form CO, and CO, in turn, is oxidized 
by both gas-phase and surface reactions to CO 2 . Actually, two cases are 
shown in Fig. I. 2. In one case, the reactor bed length is 4 cm, in the other, 

8 cm. When the bed length is cut short at 4 cm, the surface oxidation 
reactions for the hydrocarbon fuel and CO stop there. This results in a 
higher CO concentration as shown by the dotted curve. The hydrocarbon curve 
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1.0 





Figure 1.2 Detailed steady-state profiles showing the effect of catalytic 
bed length. Operating conditions: p* = 3 atm, T* = 1000°K, 

= 10 m/s, (J) = 0.15, d* = 1.8 mm and A*/A* = ^ 557 ^ 


also shows a difference after 4 cm but the difference is too small to be shown 
clearly in this figure. The gas temperature is also lower in a portion of 
the combustor for the 4 cm case. 

Fig. 1.3 gives the combustor residence time required to reach the 
emission goal vs. the residence time in the catalytic bed. The "catalytic 
bed residence time", t^, is defined by L*A*/U*^^ where is the gas 

velocity upstream of the bed and ,the gas velocity Inside the catalytic 
channel at the channel entrance is due to the area change. To 

be exact, the gas residence time in the bed should be /q 1/U^dx since the 
gas velocity in the channel, U*, is changing continuously downstream due 
to the heat release. However, since U*(x) is not known until the solution 
is achieved, the catalytic bed residence time, t*, as defined in Fig. 1.3 
is a good reference quantity and it is with this understanding when we use 
the terms of "bed residence time" and "combustor residence time". 

Referring to the curve and the calculated points corresponding to 

T*^ = 1000 K in Fig. 1.3 we see that the residence times provide an excellent 

correlation for different values of 1^*» A*/A^ and the required 

to reach emission goals. As the catalytic bed residence time, t^, increases, 

the total combustor residence time required to reach emission goals, t*^, 

first decreases sharply and then levels off. The two other curves in Fig. 1.3 

give the effects of different combustor inlet temperatures. Higher inlet 

temperature required lower combustor residence time to reach emission goals 

for a fixed catalytic bed residence time. At 1200 K inlet temperature, no 

catalytic element is needed; emission goals can be reached in an acceptable 

gas residence time for a plug flow combustor without flame stabilizer 

(t,*= o and T* = 14 seconds). This has been verified by the recent 
b EG 
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MILISECOND 


lEig. ,1.3 Minimum total combustor gas residence time for reaching emission 

goals (t*„) vs. reactor bed residence time (t*) as a function of 
£*« . 1 .. . 


combustor inlet temperature in steady-state operation, p* 

d* = 1.8 mm and T* -= 1370 K. 

af 


3 atm, 
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experimental results of Anderson (1980b). Anderson's study also suggested 
the correlation using gas residence times (1980a) . 

Transient Response 

One of the major concerns of the catalytic combustor application to 

transportation gas turbine engines is that the response time might be long. 

The computed transient results to be presented are combustor response to 

a stepwise increase of fuel flow rate (zero to a fixed value) . Its purpose 

is to guide the search for "fast-response" catalytic combustor designs. 

Some indication of the order of magnitude of the response time can be 

obtained from the definition of the characteristic time t* as defined in 

Eq (16). Small x* requires small thermal inertia, A*p*C*. If the substrate 

s s s 

density and heat capacity are fixed, small x* can be achieved with small 

solid crossectional area A*, see Fig-I»l for definition of A*. It should 

be noted that there are other parameters such as catalytic bed length and 

channel diameter which will influence the response time and they are not 

contained in x*. This will be discussed later as we see the computed results 

The combustion efficiency at the end of the combustor will be used to 

characterize the combustor response time. There are two ways to define the 

combustion efficiency. One is the so-called "carbon-balanced" efficiency, 

n [Anderson, 1975], where the combustion inefficiency is measured by the 
CB 

emission levels of CO and UHC and the energy carried away by them. The 

other efficiency is the ordinary one defined by the enthalpy rise across the 

combustor divided by the chemical energy available. If the gas specific 

heat is approximated as a constant, this efficiency, is given by 

(X* _ X* ) / (T*- - T* ) where T* is measured at the same location as the 
xn af xn 

efficiency. For a perfectly Insulated catalytic combustor in steady state. 
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and are the same. However, the two are different in transient 
CB i 

operations when there is a thermal lag in the substrate . 

For the following transient calculation, the initial surface 
temperature condition is given by T^(x,0) = 1, see Eq (21). For 1000°K 
inlet temperature, two bed lengths, 4 cm and 8 cm, are chosen for the 
transient calculation and efficiency is evaluated at their individual 
minimum combustor lengths according to Fig.1.3. Fig. 1.4 shows the combustion 
efficiencies as a function of non-dimensional time. Looking at the curves 
for (efficiency based on temperatures) , the shorter bed (4 cm) combustor 
has a faster response than that of the longer bed (8 cm) . If we take the 
time for to reach 80% as an indication of characteristic response time, 
for the 4 cm bed , t^j^ = 1.5 and for the 8 cm bed , t^^ = 2.5; namely , the 
response time of 4 cm bed combustor is only 60% of that of the 8 cm bed 
combustor. From the consideration of steady-state operation alone, one 
would choose the 8 cm catalytic bed since, as can be deduced from Fig.1.3, this produces 
the shortest required combustor length (11 cm). However, combining transient 
and steady-state considerations, the 4 cm bed combustor may be a better 
choice as its response time is 60% of that of the 8 cm bed, with the total 
combustor length only slightly longer. 

To understand why longer reactor bed results in a longer response time, 
surface and gas temperature distributions are plotted in Fig, 1.5 as a 
function of time. From the substrate temperature distribution we see that 
during start-up transient, the solid temperature is higher in the flow 
entrance region as a result of higher mass transport and surface reaction 
rates. Comparing the solid and gas temperature histories for the case 
when catalytic bed length is 8 cm, we see that at t = 0.5, the solid is 
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Fig. 1.4 Combustion efficiencies as a function of time during start-up 
transient, p* = 3 atm, (f) = 0.15 and A*/A* = .667. 




X* ( CM ) 


X* ( CM ) 


Fig. 1.5 


Gas and substrate temperature history during transient start-up 
for two bed lengths. Same operating and design parameters as 
those in Fig. 1.4. 





transfering heat to the gas in the upstream region in the bed. At 
X* = 3.5 cm, the gas temperature becomes higher than the substrate temperature, 
heat is then transferred from the gas to the solid in the down-stream portion 
of the reactor channel. If the heat loss rate is greater than heat generation 
rate due to gas phase reactions, then the gas temperature will drop as shown 
in Fig.1.5. If we cut the reactor bed short at 4 cm, a portion of this heat loss 
can be eliminated and the downstream gas temperature (x* > 4 cm) will be 
higher as indicated by the dotted curves in Fig.1.5. This explains why a 
shorter reactor bed can have a faster response. 

More computation of the response time dependence on the catalytic bed 
length, gas velocity and combustor inlet temperature is shown in Fig. 1.6. 

As will be shown later, the nondimensional transient response also correlates 
well with the gas residence times in the reactor bed and the after-bed space. 
Therefore the results will be presented using residence times rather than 
length and velocity. In Fig. 1. 6, the adiabatic flame temperature, pressure, 
channel diameter, open area fraction and are held constant. The 

solid line in Fig. 1. 6 gives the response time if the reactor-bed-gas-residence 
time is the minimum to reach emission goals in steady-state operations. When 
the combustor inlet temperature is increased, the minimum required reactor- 
bed-residence time is decreased and the combustor response time is shortened 
drastically. If a more conservative design is used, i.e. longer bed residence 
time than is required, the combustor response time is increased for a 
fixed inlet temperature. The increase is very slight in the beginning when 
the bed residence time is close to the optimum value but becomes larger as 
the bed residence time increase further. This implies that a conservative 
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Nondimensional response time (n,p goes to 80,^) as a function of 
combustor inlet temperature and reactor bed residence time. 
L*/U*e£ = 15 msec, T*^ = 1370 K, A*/A* = .667, d* = 1.8 mm and 

p* = 3 atm. 


Fig. I. 6 




design with reactor bed length somewhat longer that that of the optimal one 

will not greatly increase the combustor response time. 

Fig, I. 7 gives the nondimensional response time as a function of channel 

hydraulic diameter and flame temperature. For all the calculations in Fig.1.7 

the residence time in the catalytic bed, t*, is held constant at 1 milisecond 

b 

and the residence time in the after-bed space, t* , is 13,5 miliseconds. The 

^ as’ 

computation for 1370 K and 1590 K includes different combinations of L*, L* , 

3.S 

U*ef and A*/A* and it can be seen from Fig-I»7 that the nondimensional response 
time correlates well with the residence times. It can, therefove, he oonotuded 
that the primary influence of the catalytic hed and after-hed lengths^ the 
approaching gas velocity and the open area percentage is contained in the gas 
residence times in the catalytic hed and in the after-hed space for both steady 
and transient operations. 

Fig. 1.7 also shows that for a given flame temperature, the nondimensional 
response time increases with decreasing channel diameter d*. This, however, 
does not necessarily mean that the dimensional response time increases with 
decreasing d*. If the open area percentage (or equivalently A*/A*) is 
fixed, decreasing d* decreases which is used in the reference time scale 
as defined in Eq (16) and the dimensional response time may actually decrease. 
Table 1.3 gives selected dimensional response times for ceramic and BCanthal 
metal alloy substrates using the computed results from Fig.1.7. it can be seen 
that in general the dimensional response time decreases with decreasing 
channel diameter but the amount of decrease vanishes when the channel diameter 
becomes too small. In special cases such as the 1590 K flame temperature 
case, the response time actually reaches a minimum and increases with further 
diameter decrease. 
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Fig. 1.7 Computed nondimensional catalytic combustor response time 

correlation with gas residence times for various channel 

diametefs and flame temperatures, p* = 3 atm, tf = 1 msec, 

t* = 13.5 msec, T* = 1000 K. 
as in 
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Table 1-3 Computed Catalytic Combustor Response Time Using Ceramic 

and Kanthal Metal Substrates. 

t* = 1 msec, t* =13.5 msec, p* = 3 atm, T* = 1000 K 
b as ' in 

t* (sec) t* (sec) 

r r 

(K) A*/A* d* (mm) ceramic Kanthal alloy 

3.x S 


1370 

10 

1.8 

.848 

1.52 

1370 

10 

1.0 

.451 

.807 

1468 

10 

2.5 

.909 

1.63 

1468 

10 

1.8 

.533 

.956 

1468 

10 

1.0 

.454 

.813 

1590 

10 

2.5 

.527 

.944 

1590 

10 

1.8 

.377 

.675 

1590 

10 

1.0 

.448 

.802 

1590 

5 

1.8 

.754 

1.35 

1590 

2 

1.8 

1.88 

3.37 
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If the gas residence times and channel diameter are held constant, 
varying the open area percentage (or A*/A*) changes the dimensional response 
time by a factor of A*/A*. This is the direct consequence of the correlation 
between nondimens ional response time and residence times as just discussed. 
Large open area percentage (or thin substrate), therefore, favors fast response. 

Fig 1.7 indicates higher flame temperature shortens the response time 
for large diameter channel but flame temperature makes little difference in 
response time for very small diameter channels as all curves converge. 
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CHAPTER II STEADY STATE COMPUTATION AND COMPARISON WITH EXPERTMF,NTS 


Three-step Semi-global gas-phase Chemical Reactions 

Gas phase reaction studies between hydrocarbons and air in lean mixtures 
have shown (Dryer and Classman, 1978) that the hydrocarbon fuel first decom- 
poses into a number of lower hydrocarbons with the simulataneous formation 
of some water vapor. The combination of these "pyrolysis" steps results in 
an almost isothermal region after which CO is formed and temperature begins 
to rise. The carbon monoxide is then oxidized and the reactions are completed. 
These kinetic studies are very useful for the development of semi-global 
kinetic schemes. 

Two recent catalytic combustion models, by Bracco et al (1980) and T'ien 
(1980), have utilized the gas phase semi-global modeling concept. The study 
by T'ien was actually for transient response but a quasisteady gas phase 
was assumed. The gas phase reactions consist of two semi-global steps; the 
hydrocarbon reacts with oxygen to form CO and H 2 O and CO then reacts with 
oxygen to form CO 2 . Oii the surface, both hydrocarbon and CO react with oxygen 
to form CO 2 and or H 2 O. The steady-state prediction of that model agrees 
qualitatively with the experimental trend of CO and hydrocarbon emissions, but 
a quantitative comparison has not been made. From a pure gas phase reaction 
point of view, two-step semi-global scheme predicts too high a level of CO 
in the initial pyrolysis region even though the CO peak level and the decay 
curve can be made to agree with the experimental data by kinetic parameter 
adjustment. This early CO appearance problem can be alleviated to a certain 
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extent by Increasing the number of appropriate semi-global gas reaction 
steps. In the work by Bracco et al (1980), for example, a three-step 
reaction scheme was used for propane oxidation. Propane (C_Hq) was assumed 

J o 

to form C 2 H^ first, ^211^ to CO next and then CO to CO 2 . 

For most practical fuels, the pyrolysis intermediates are numerous, and 
a dominant intermediate species may or may not exist. Since in high inlet 
temperature operations (T*^ > 800*K) , the catalytic combustor length required 
is limited by the emission level of CO, not hydrocarbons (Anderson, 1980a), 
it seems reasonable to develop a scheme which gives the correct CO character- 
istics without being tied down by the uncertain nature of the hydrocarbon 
intermediates . 

In this chapter, the steady-state computation will be based on a three- 
step gas-phase kinetic scheme. The semi-global reactions are; 


(A’) C H 

^ n m n m 


d t* 


= - cy^i[c H ]'2T"’3exp(-E*/RT*) 




u 


n m 


( 22 ) 



“4exp(-E*/RT*) 

d t* ' " ( 23 ) 


( 
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(B) CO + I 0, -> CO 

C 2 2 


mi 

d t* 


C*p*^[C0]^2[02]^3[H20]^5T*^ X 

exp (-E*/RT*) 


( 24 ) 


In reaction (A'), Eq . (22), the original hydrocarbon fuel forms the 
"pyrolyzed hydrocarbon", in an isothermal reaction. The "pyrolyzed 

hydrocarbon" then reacts with oxygen to form CO and H^O in reaction (A") . 

In reaction (B) , CO is oxidized to CO 2 . The introduction of reaction (A') 
serves to decrease the amount of early CO appearance as found using the 

(p) 

two-step model. The identity of is of no great consequence for the 

purpose of this modeling work, if the mass diffusion coefficient is properly 

(p) 

evaluated. For convenience, in the present calculation we take C H to 

’ n m 

have the same n and m as the original fuel. 

Corresponding to the reactions given by Eqs. (22-24), the species 
equations are : 


Original Hydrocarbon: 

“dx VhC "VhC,s^ " 


pu 


w 


(25) 
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"pyrolyzed" Hydrocarbon: 
pu 

Carbon monoxide: 



Wq-w 


1 


PU ^ ^ , W,C^ - W2 

Oxygen : 
where 

VhC " ^HC 


The reaction rates are given by 
Wo = BoP’'>p''n''3 YhI «p(-Eo/T) 
w,=B,p“lp(“2 + “3 )t»u •^“2 exp(-Ei/T) 

W 2 = B 2 p«lp(«2+63*e5)T^ y ^2 yn Y^es exp(-E 2 /T) 
and 

j Y^(,(o) 


( 26 ) 

( 27 ) 

( 28 ) 

( 29 ) 

( 30 ) 

( 31 ) 

( 32 ) 
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Three global catalytic surface reactions are assumed to occur on the 
surface. They are: 

(C) + (n4^) O 2 -> nC02 + | H 2 O 

with the corresponding surface rate J* given by 
[C^HJ 3 exp (-Ef/RTp 


(c") ^ >^*^02 + 2 H 2 O (34) 

with j; = cUc„H^P>]j exp (-Ej/RTp 

(D) CO + J O 2 CO 2 (35) 

with J* = C*[C0]3 exp (-E*/RTp 


Other equations are similar to those in Chapter I. The steady solution 


is obtained from the transient combustion code as discussed in Chapter I. 
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II. 2 Computed Results 


The design parameters of a catalytic combustor include the catalytic 
bed length, the after-bed reaction distance, the channel hydraulic diameter, 
the open area percentage, the catalyst variables (catalyst t3q>e, loading 
density, washcoat type) and staging arrangements (graded cell, for example). 

The operating parameters include inlet temperature and pressure, the 
approaching (reference) mixture velocity, the fuel/air equivalence ratio 
or the adiabatic flame temperature and fuel types. An exhaustive parametric 
study will not be practical for this report. Instead, selected computations 
are presented. Since emphasis is placed on automotive gas turbine application, 
we are most interested in higher inlet temperature, above atmospheric pressure 
and high reference velocity operations with thin substrate and large open area 
percentage design. The total combustor length (bed plus after-bed) required 
to meet the emission goal (13.6 g CO/kg fuel, 1.64 g THC/kg fuel) as a func- 
tion of these parameters is the main output of the model calculation. 

The chosen fuel is No. 2 Diesel which was used by Anderson (1977a, 1978, 
1980a). Its hydrogen-carbon atom ratio is 1.8 with 27.55% aromatics. The 
molecular weight is taken to be 138 for the purpose of estimating the diffusion 
coefficient which is assumed to vary as one over the square root of molecular 
weight. Other properties can be found in Anderson (1977a). The catalytic 
combustor used for reference calculation is the one by Matthey Bishop which 
was also described by Anderson (1980a) (MBI-2.5 and MBI-5.4). The values of 
kinetic parameters are listed in Table II. 1. Expressions in Reactions A" to D 
are taken from Edelman and Fortune (1969), Dryer and Classman (1973), Marteney 
(1979) and Kuo and Morgan (1971) respectively, all with adjusted pre-exponential 
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Table 11,1 Chemical Kinetic Constants 
(for Chapter II) 


Reaction 

(A') 

Reaction (A") 

Reaction (B) 

Reactions (C')(C") 

Reaction (D) 

C* = 1.6 X 
o 

10 ^ 

C* = 1.8 X 10^ 

^1 14 

C* p* = .54 xlO 

C* = 1.5 X 10^ 

C* = 10^ 
4 




(at 3 atm) 



Yl = 0 


= 0.3 


E* = 10 
(kcal/mole) 

E* = 17.8 
(kcal/mole) 

Y 2 = 0.5 


= 0.5 

62 = 1 



II 

o 


“3 = ^ 

8 ^ = 0.25 





“4 = ^ 

= 0 



E* = 34 


E* = 24 

65 = 0.5 



(kcal/tnole 

) 

(Iccal/mole) 





E* = 40 


(kcal/mole) 



factors. The calculated results are not very sensitive to catalytic kinetic 
parameters indicating mass transfer dominates the surface reactions. The 
choice of gas-phase hydrocarbon oxidation kinetic constants deserves some 
explanation, since kinetic rates for No. 2 Diesel fuel are hard to find. 

Ignition delay times as a function of initial temperature, however, do 
exist for long chain and cyclic hydrocarbons (Dryer and Classman, 1978, 

Edelman, 1978), which are the ingredients of the Diesel fuel. The ignition 
delay time of No. 2 Diesel is assumed to be the arithmetic mean of its 
ingredients. Model calculations are then performed for pure gas phase reac- 
tions and the ratio of calculated hydrocarbons disappearance times as a 
function of initial temperature is taken to be proportional to the ratio of 
ignition delay times as indicated in Fig. II. 1. This provides some measure 
for the correct temperature dependence in the gas kinetic parameters for the 
diesel fuel. The use of "hydrocarbon disappearance time" in this comparison 
reflects the difficulty of identifying ignition time in the theoretical 
study. The pre -exponential factor in reaction (B) is adjusted to fit the 
CO decay data of Anderson (1980a) . 

The model calculations for the catalytic combustors will be presented next. 

Effects of Catalytic Bed Length, Downstream Reaction Distance, 

Reference Velocity and Adiabatic Flame Temperature 

Figure II. 2 shows the detailed profiles of two cases; one with a 5-cm bed 
and the other with a 2. 5-cm bed. If we look at the 5-cm case first, we see 
when mixture enters into the reactor channel, the fuel is quickly consumed 
on the surface by catalytic reaction, the surface temperature starts from a 
low value (T^ = 1.106), due to the fast heat and mass transfer rates at the 
entrance region, and rises sharply away from the channel entrance. The gas 
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Fig. II. 1 Comparison of temperature dependence for the computed disappearance 
time (nondimensional ) of Diesel No. 2 fuel with Ignition Delay 
Data (Edelman, 1978). Reference time is taken at 1200 K. 
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Fig- II. 2 Computed temperature and species profiles for two similar cata- 
lytic combustors with different bed lengths^ A*/A| = 0.93 and 
d* = 1 .382 mm. 
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temperature begins to rise due to the heat transfer from the substrate and 
when T reaches 1.15, gas-phase reaction is initiated and the hydrocarbon fuel 
is converted to the "pyrolyzed hydrocarbon,”, which is equal to the difference 
between (total hydrocarbon) and (original hydrocarbon) . At the 

same time, CO, which is produced only in the gas-phase, begins to rise. The 
carbon monoxide in turn is consumed by both gas-phase and surface reactions 
to form CO 2 until the bed exit is reached. After the exit, CO is consumed 
by gas-phase reaction only so a slope discontinuity of the CO curve occurs 
at X* = 5 cm. Similar discontinuity can be found in temperature and hydro- 
carbons, although in the latter case, it is not very apparent. If at the 
bed exit the gas temperature is high enough, the gas-phase reaction will 
continue at a high rate and adiabatic flame temperature is reached at x* = 12 cm. 
As shown in Fig. II. 2, a CO peak occurs when its production and oxidation rates 
become equal. 

If the catalytic bed is cut short at 2.5 cm. as shown in Fig. 11.2, the 
total hydrocarbon ( ^ oxidation rate decreases suddenly at the bed 
exit (2.5 cm) due to the absence of surface reaction. Since the gas 
temperature is not yet high enough, the gas-phase reaction downstream 
proceeds slowyly as shown by the slower Increase of gas temperature. Because 
a larger fraction of the hydrocarbon fuel has to be consumed in the gas- 
phase reaction, higher CO peak occurs and the total combustor length 
required to reach the emission goal is accordingly longer than the case 
with longer catalytic bed. 

It was shown in Chapter I that the effects of catalytic bed length, 
after -bed distance, gas velocity and open area percentage could be correl- 
ated using reference gas residence times in the reactor bed and in the after- 
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bed space. For a fixed open area fraction, this correlation is simplified 

and can be expressed by the ratios of lengths over the reference velocity 

as given in Fig. II. 3. For a fixed "gas residence time" in the bed (L*/U*^j) , 

there is a corresponding minimum combustor residence time for 

EG ref 

emissions goals to be satisfied. Although the residence time correlation is 
natural for gas-phase reactions, it is less obvious for surface reactions. 

The validity of this type of correlation is independent of the gas kinetic 
schemes used as illustrated here and in Chapter I where two-step kinetics and 
different sets of catalytic parameters were employed. 

Fig. II. 3 shows that as the catalytic bed length increases, the required 
combustor length decreases. This may lead one to choose a design which is 
all reactor bed without the downstream gas reaction space. However, as 
discussed in Chapter I, in many catalytic combustor applications (e.g., 
transportation engines), a fast combustor response is required. Longer 
reactor bed results in a longer response time due to the larger substrate 
thermal inertia. Therfore, certain compromises may have to be made between 
steady and nonsteady considerations. Take Fig. II. 3, for example, one may 
want to choose a bed residence time of 4 milliseconds instead of 6 milli- 
seconds. This increases the total combustor length by 15% but reduces the 
response time by 30%. As the inlet temperature becomes higher, as preferred 
by the automotive gas turbine design, the required reactor bed will be 
drastically reduced, leaving most reactions to occur in the gas-phase. A 
shorter reactor bed also decreases the pressure loss. 

The importance of downstream gas-phase reactions has been recognized 
by Anderson as demonstrated in much of his work (Anderson, 1977b, 1980a, T'ien 
and Anderson, 1979). Figure II. 4 compares the model calculation with his re- 
cent data (Anderson^ 1980a) on CO emission index as a function of downstream 
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Fig . II. 3 Computed "residence time" correlation: total residence time 


in combustor (bed plus after-bed) to reach emission goals vs 
residence time in catalytic beds. 
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Fig . II. 4 CO emission index vs. adiabatic flame temperature: 
of model calculation with Anderson's data (I980a) . 


comparison 
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reaction distance and adiabatic flame temperature. Figure II. 5 compares the 
minimum operating flame temperature vs. the downstream reaction distance for 
two bed lengths and two velocities. The comparison between the model and 
data is reasonably good within the range of parameters varied. 

The appearance of CO peak as indicated in Fig. II. 2 has been reported 
in previous experimental work (T'ien and Anderson, 1979, Bracco et al, 1980). 

Effects of Inlet Temperature and Pressure 

Figure II. 6 shows the computed combustion efficiency as a function of 
adiabatic flame temperature for three inlet temperatures and two pressures 
at three different downstream locations. For a fixed pressure and adiabatic 
flame temperature, higher inlet temperature results in a higher efficiency 
at all three downstream positions. This is consistent with the experimental 
trend found by Anderson (1978). With higher inlet temperature, the gas- 
phase reaction is initiated earlier which results in a higher efficiency at 
a fixed location. However, as more residence time is given to the mixture, 
the difference in efficiency diminishes as can be seen from curves at differ- 
ent locations in Figure II. 6. 

From Fi& II. 6, for inlet temperature of 1000 K, the combustion efficiency 
at 3 atm is always higher than that at 6 atm. This is because of the faster 
mass diffusion rate at lower pressure. For 1100 K inlet temperature, however, 
a portion of the 6 atm efficiency curve is higher than that of the 3 

atm case as can be seen in Fig. 11.6(b) and (c) . This reversal phenomenon 
was first reported by Anderson (1978) who suggested that this was due to the 
competition between gas-phase and surface reactions. This model calculation 
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Fig. II. 5 Minimum operating flame temperature (to reach emission goals) 
vs. downstream reaction distance: comparison of model calcu- 

lation with Anderson's data (I980a) . 
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Fig. II. 6 Combustion efficiency vs. adiabatic flame temperature for 

different inlet temperature and pressure at three downstream 
locations. 
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supports this interpretation. High temperature and pressure support 
fast gas-phase reactions, but high pressure also retards surface reaction 
through slower mass diffusion rate. Taking the case of 1100 K inlet tempera- 
ture and 1320 K adiabatic flame temperature for example: at the bed exit, 

the 6 atm case has a lower efficiency than the 3 atm case because inside the 
reactor channel, surface reaction dominates. At downstream distance of 5 cm, 
the 6 atm case still has lower efficiency. But at 10 cm downstream, the two 
efficiencies become equal and at 15 cm downstream 6 atm has a higher effi- 
ciency. Since only gas-phase reaction exists in the after -bed space, the 
efficiency reversal occurs when the gas-phase reactions overpower the upstream 
influence of surface reactions. 
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SUMMARY 


A transient model of monolith catalytic combustor is presented in this 
work. The model assumes a quasisteady gas phase and a thermally-thin solid 
with the substrate thermal inertia responsible for the combustor response 
time delay. In the gas-phase treatment, several quasi-global chemical 
reactions are assumed which are capable of describing CO and unburnt hydro- 
carbon emissions in fuel-lean operations. By neglecting heat conduction 
along the flow direction in the gas and the solid in high speed approaching 
flows, the system of differential equations describing the combustor transient 
is simplified and Runge-Kutta integration scheme is applicable. The resulting 
computation scheme is highly efficient in computational time and is suitable 
for parametric calculation for both steady states and transients. 

In the steady-state computation presented, the influence of selected 
operating and design parameters on the minimum combustor length is studied. 

Special attention is given to the effect of after-bed gas phase reaction 
space. Comparison with the data by Anderson indicates that the model is 
able to describe all the salient features found in the experiments including 
the appearance of CO peak and the efficiency reversal phenomena at high 
inlet temperature and pressure. Quantitative agreement with CO emission 
data is possible if a certain degree of rate constants adjustment is made. 

The model computation also suggests that for determining the condition of 
achieving the steady-state emission goals, the gas residence times in the 
catalytic bed and in the after-bed space are the similarity parameters 
absorbing the influence of catalytic bed and after-bed lengths, the gas velocity 
and the open area fraction of the catalytic element. Likewise, the nondimens ional 
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response time and the gas residence times are the proper similarity parameters 
in transient operations. A consequence of this is that the dimensional 
response time is proportional to the substrate crossectional area A* (all 
other parameters fixed) and so a thin substrate is necessary for fast response. 

The computed time history provides an understanding of how a catalytic 
combustor responds to an upstream condition change. The leading edge of 
the substrate is the first to respond to an upstream variation, say, a change 
of fuel flow rate. The leading edge is closer to the source of variation 
and the mass transfer rate is higher in the flow entrance region of the gas 
channel. As a result, the surface fuel concentration is perturbed first 
which, through the surface reaction, produces a surface temperature change. 

Heat transfer across the flow channel then causes the gas temperature to vary 
and this disturbance propagates downstream by convection. Since, the down- 
stream part of the substrate always responds to the upstream substrate 
temperature variation, it is not difficult to see from the above description 
that a combustor with a longer catalytic bed residence time will take a 
longer time to reach a new equilibrium state. 

The model calculation also shows that higher combustor inlet temperature 
and adiabatic flame temperature shorten the response time. Smaller channel 
hydraulic diameter usually results in a faster response but there are exceptions . 
Since the requirements for faster transient response do not always coincide 
with the optimal steady-state design conditions, for applications where 
transient characteristics are important, both steady and transient studies 
have to be made simultaneously so that a good compromising design can be 
found . 
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Nomenclature 


A* 




A* 


crossectional area of one gas channel 

substrate cross-sectional area associated with one gas channel, 
see Fig. 1 


= A* + A*, total crossectional area of one cell unit 


k?) Y^^c(o) 


W, 


0-^2) 

HC 


L* 


1 PA t\A r\A ^4 

tttttt ^(0.0^ to.o) ^to.o) 


u*(0.0) p*(0,0) Yjj^(O.O) 


„ y (1-02) W 

* HC 02 


u*(0.0) p*To.oT €7(0,^ 2 p*(o.o) p*(o,o) t*(o,o) 

X Y®3,„ 0, „ (1-62) „ -63 „ -65 

X YQ2iO,OJ Wq2 


= C* d* / a* Nu^ 


B, = C* d* / a* W-- Nu 
4ft CO •• 


" '^CO ^ ^HC 


C* 

s 

C* 

P 

d* 


d* 

s 


specific heat of substrate 
gas specific heat 

catalytic monolith combustor channel equivalent diameter 
half thickness of tbo substrate 
activation energy 


mass consumption rate of hydrocarbons per unit surface area 
through reaction (C) 


Di 


mass consumption rate of CO per unit surface area through 
reaction (D) 


(Lep 


2/3 


'H 


4 Nu (o*(0,0)/d* )(L*/u*(0,0)(k*/k*(0,0) 
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k* 

L* 



L* 

c 


L* 

EG 


Le, 

ID 

n 

Nu 


Nu 


P 

'll 



S* 

T 

T 

s 


heat conductivity of the gas 
catalytic bed length 

length of after-bed downstream gas phase reaction space 


total combustor length (catalytic reactor length plus 
after-bed downstream gas reaction distance) 


total combustor length (catalytic bed plus after-bed space) 
required to meet emission goals 

Lewis number for hydrocarbon gas in air 

Lewis number for CO in air 

number of hydrogen atom in C H 

n m 


number of carbon atom in C H 

n m 


Nusselt number at x 

Nusselt number for fully developed flow 
p*/p*(0,0), nondimensional pressure 

q*/C*T*(0,0) , nondimensional heat of combustion per unit mass 
of C H in reaction A 


n m 


q*/C*T* (0 ,0) , nondimensional heat of combustion per unit mass 
of CO in reaction B 


q*/C*T* (0 , 0) , nondimensional heat of combustion per unit mass 

of C H in reaction C 
n m 


q*/C*T*(0,0), nondimensional heat of combustion per unit mass 
of CO in reaction D, same as q^ 
see Eq. (17) 

circumferential length of one channel crossection (* ird*) 
T*/T*(0,0), nondimensional gas temperature 
T^*/T*(0,0), nondimensional substrate temperature 
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adiabatic flame temperature 


T*, 

af 


T* 

in 




t* 


as 


combustor inlet temperature 
m nondimenslonal time 

= (L*/U*^P (A*/A*), reference gas residence time in the catalytic 

combustor bed 

= L* /U* reference gas residence time in the after-bed space 


t*^ characteristic transient response time (^^->-80%) 


t* 

EG 


ref 


W, 


w. 


Y . 
1 


'^CO 

^02 

%c 

dHC 

^THC 


Minimum total combustor gas residence time (t* + t* ) needed to 

D as 

reach emission goals 

u*/u*(0,0), nondimenslonal flow velocity 

reference gas velocity, neaaured at upstream of catalytic bed 
molecular weight of species 1 
chemical reaction rate 

x*/L*, nondimenslonal axial distance measured from catalytic bed 
entrance to downstream 

mass fraction of species 1 






^HC/ 


y +L(p) 

^HC ^HC 
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a* thermal diffusivity of gas 

rate exponent, see Eqs (1) and (23) 

6 rate constant, see Eq. (22) 

rate exponent, see Eqs (2) and (24) 
rate exponent, see Eq (22) 

p* gas density 

p* substrate density 

s 

T* characteristic substrate heat-up time, see Eq. (16) 

Ti__ carbon-balanced efficiency 

LJd 

combustion efficiency based on temperature difference 
1 concentration, g-mole/c.c. 


Subscript 

s surface 

Superscript 

^ dimensional quantity 
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APPENDIX I 


LISTING OF COMPUTER PROGRAMS FOR CHAPTER I 


The numerical computations were performed using an interactive VAX - 
11/780 computer. FORTRAN language was used in writing the main 
program (e.g. IGNITION. FOR; 17) and the data program (e.g. BLOCK. 

FOR; 31). They were then compiled to form the objective programmes 
(e.g. IGNITION. OBJ; ** and BLOCK. OBJ; **) . The two objective 
programmes were then linked together to form an executable program 
(e.g. IGNITION. EXE; 10). In performing a calculation, only the 
executable is needed . 

The BLOCK program stores most of the data which are not likely 
to be changed from one calculation to the other, while the most 
frequently varied input data will be read in in each computer run. 

If the data stored in the BLOCK program are to be changed, the FORTRAN 
BLOCK program has to be modified (edited) , re-compiled and linked 
with the main program to form a new executable program. 

In the following, the FORTRAN programmes (IGNITION. FOR; 17 
and BLOCK. FOR; 31) are listed together with a sample calculation 
showing input and output data. 

The IGNITION program is presently set up to perform the start-up 
transient. With minor modification, it can read in any prescribed 
initial condition and to calculate the corresponding unsteady response. 
Given enough computational time, the IGNITION program will also 
result in the steady-state profiles. 

More detailed instrictions (for IGNITION. FOR; 17): 

1. The number of grid points used should be less than 300. 

If more than 300 points are needed, line 600 has to be 
changed . 

2. Read-in data (lines 10900 - 12300): 

XPHI = (J)= fuel-air equivalence ratio 
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TI = T (0,0) = upstream, initial temperature (K) 

XPT = P(0,0) ^ upstream, initial pressure (Atm.) 

UIUS = “ upstream (outside oT the bed) Initial gas 

velocity (m/s) 

XLS = reference length (m) , suggested to use 0.1 m (10 cm) 

DS = d* = channel hydraulic diameter (m) 

AOAS = A*/A*^ = open-to-close area ratio 

Djj = nondimens ional step size in X (space) . It is 

nondimensionalized by XLS. So if XLS = 0.1 m and 
DX = 0.01, dimensional step size is 0.001 m or 
1 mm. 

DT = nondimensional time step size 

TMAX = maximum nondimensional time allowed for computation 
IPRINT = time interval for printing out transient profiles. 

For example, if DT = .01 and IPRINT = 25, 

Output will be printed out every t = 0.25. 

N = number of grid points (in X direction) in the 

catalytic bed. The nondimensional catalytic bed 
length is given by (N-1) x (DX) . For example, 
if N = 51, DX = .01, then the nondimensional bed 
length is equal to 0,5. If XLS = .1 m, then the 
dimensional bed length is .05 m (5 cm). 

NAS = number of grid points (in X-direction) in the 
downstream after-b .d space. The length of the 
downstream distance can be calculated in a similar 
manner as in the above. 

Cl, C2, C3, C4 = C^, C^, C^, (See Table 1.2) 

IT = interval in X for printout. For example, if DX=0.01 

IT = 5, the printout will be at X= 0, 0.05, 0.1, 0.15, . 
etc. 

3. Those data which are not read in will be defined by the BLOCK 
program. 
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4. Some Symbol Definitions: 

PR = Prandtl number of air 

REY = Reynolds number based on channel diameter and 
upstream velocity 

ES,S = E^*, E 2 *, E^*, E^* (See Table 1.2) 

E,S = E^, E 2 J , E^ 

C,S = C^*, C 2 *, C^*, C^* (See Table 1.2) 

GC = gas constant (joule/kg K) 

XLEWl = Lewis number for C H in mixture 

n m 

XLEW2 = Lewis number for CO in mixture 

3 

ROI = initial, upstream (I.U.) gas density (kg/m ) 

XKI = I.U. gas heat conductivity (cal/m sec K) 

CP = specific heat of gas at constant pressure (cal/kg K) 
ALPHA = I.U. gas diffusivity (m /sec) 

™ - Vo, 

TMAX = Maximum time for computation 

5. Output data (lines 40300 to 42200) 

The following quantities are printed out as a function 
of axial distance and time: 

Nondimens ional gas and surface temperatures, Emission 
Indices (xlOOO) of hydrocarbon fuel and CO (both in the gas 
phase and on the surface) and carbon-balanced and thermal 
efficiencies . 


TS = 

T 

s 

YES = 

5^HC,s 

YCS = 

f^C0,s 

YF =] 
c 

^HC 

YC =t 

4 

^CO 
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IGNITION. for; n 


12-APR-1981 23:29:54.94 


PAGE. 1 


100 
200 
300 
■•00 
500 
bOO 
700 
100 
900 
1000 
1100 
1200 
1300 
1400 
1500 
1600 
1700 
1800 
1900 
2000 
2100 
2200 
2300 
2400 
2500 
2600 
2700 
2750 
2800 
2900 
3000 
3100 
3200 
3500 
3600 
3700 
3800 
3900 
4000 
4100 
4200 
4300 
4400 
4500 
4600 
4700 
4100 
4900 
5000 
5100 
5200 
5300 
5400 
5500 
5600 
5650 
5700 
5800 
5900 
6000 
6100 
4200 
4300 
6400 
6500 
4600 
4700 
4100 
6900 
7000 
7100 
7200 
7300 
7400 
7500 
7600 
7700 
7100 
7900 
1000 
• 100 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


7 

• 

9 

10 

11 

12 

13 

14 

15 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

31 

32 

33 

34 

35 

36 

37 

41 

42 

43 

44 

45 

46 

47 
41 
49 

71 

72 

73 

74 

401 

402 

403 

404 

405 

406 

407 
401 

409 

410 


♦J*J4»»CAIALniC^CCJEIJST0R^STAPT-UP^TRANSIENT44»»44»»4*44»'»4» 

PARAMETER NX=300 

COMMON /ACTI/ESl ,ES2,ES3,ES4 

COMMON /MOLEO/ WHC, WCQ, *02, Wh20 

COMMON /THI/XN, XM, CIS, 025 

COMMON /TH2/ CP, GC, XPH 

COMMON /UPST/ VCl, YCI 

COMMON /1N1/TS,1ES,YCS 

DIMENSION YF(NX), TCNX), TS(NX), YO(NX), YFS(NX), KC(NX), U(NX) 
DIMENSION XKU4), XKTt4l, XKTSC4), XKCC4), YC(Nx5, YCSIN)() 
DIMENSION XJH(NX), XJDKNX), XJD2(NX), YH2C(NX) ,XP(NX) 
dimension EFFCNXl, RNN ( N X ) , EFF T ( NX ) 

FORMAT (4F20.10) 

FORMAT (14) 

FORMAT (IHO.'SCLID TEMPERATURE (TS)V) 

FORMAT (6F12.6) 

FORMAT (IHO, 'SURFACE HC MASS FRACTION (YFS)'/) 

FORMAT (//) 

format (IHO, 'GAS TEMPERATURE (T)'/) 

FORMAT (IHO.'CNHM MASS FRACTION (YF)'/) 

FORMAT! IHO, ' Y 0 ( UPST, INI ) =',F10.6,5X,'YF(UPST,INI)s',F10.7 
1,5X,'E0U1 RATIO ='F15.9//) 

FORMAT (1HO,’01=',F15.7,5X,'C2=', F15.7,5X,'DX=',F15.7,5X,'CT=' 
lFl5.7//) 

FORMAT (IHO, 'N=' ,110//) 

FORMAT ( IHO, ' INITIAL CONDITIONS (TS)',/) 

FORMAT (1H0,'TMAX = ',F15.8) 

format (IHO, 'NO, OF GRID POINT =’,2110//) 

FORMAT (IHO, 'TIMES ' ,F 15. I,5X, 'MOUNTS' ,110/) 

FORMAT (///) 

FORMAT (10E12.5) 

format (IHO, 'INITIAL CONDITIONS (YFS)',/) 

FORMAT (lHO,'XJh'/) 

FORMAT dhO,' CC MASS FRACTION (YC)'/) 

FORMAT (IHO, 'SURFACE CO MASS FRACTION (YCS)'/) 
format (IHO, 'INITIAL CONDITIONS (YCS)',/) 

FORMAT (IHO, ' ES,S s ',6F12.b/) 
format (IHU, ' C,S s ',4E20.10/) 

FORMAT (IHO, ' A,S s ' ,6F12.6/) 
format (IhO, ' PE,S = ',6Fl2.b/) 

FORMAT (IHO, ' YCO =' ,/) 

format (IHO,' YCO(SURFACE) =’,/) 

format (IHO,' L/U s', F15.7/) 

FORMAT (IHO,' PIS, B2S, s', 4E15.8/) 

FORMAT (IHO.'Bl, B2, 63, B4 s',4£15.8/) 

FORMAT (4E20,8) 

FORMAT (IHO,^ E,5 s',6F12.6/) 

FORMAT (1H0,'GC, XLEWl, XLEW2, XLS, PS =',6F12.6/) 

FORMAT (1H0,'RCI, XKl, CP, ALPHA s » , 3F 1 2 . 6 , E 1 2 . 5 / ) 


FORMAT (IHO, 'APPROACH VEL = ',FiO,4/) 

FORMATdHO, 'CHANNEL VEL d.U.) * ’,F10.5/) 

FORMAT dHO,'EFFlCIENCY(CARBCN-BALANCE)',/) 

FORMAT (IHO, 'EFFICIENCY (TEMPERATURE) ' ,/) 

FORMAT dhO,' XFI s', F20.5/) 

FORMAT (1H0,'TI, PR, REY, AD2UL, ADAS s',5F12.5/) 
FORMAT dHO,' ES,Ss'4E20.8/) 

FORMAT dH0,'4444*»444»4*»4*4»4»**A»»*»»*»A' ,/) 
FORMAT (IHO, 'GAS TEMPERATURE IS NEGATIVE') 

FORMAT ClHO,' CNHM CONCENTRATION 15 NEGATIVE') 
FORMAT dHO,' CO CONCENTRATION IS NEGATIVE') 
FORMAT dHO,' YFF IS NEGATIVE') 

FORMAT (IHO, 'ADIABATIC FLAME TEMP s' F10.6/) 
FORMAT (1H0,'AD1AB FLAME TEMP(K) e',F15.5/) 


Sb««$b«S**b4«*»*»*4*«**4*«*4**«4«4**«bbt*« 

»«b«s*«s*b«»XNPUT PARAMETtBS»4*»*»4*»»»»»' 
•»4»b*»»»*CHEMICAL PARAMETERS************ 

CHEMICAL REACTICN 1, CNHM ♦ C2 s CO ♦ H2C (GAS PHASE) 
CHEMICAL REACTICN 2, CO + 02 s C02 (GAS PHASE) 

CHEMICAL REACTICN 3, CNHM ♦ 02 s C02 4 H2C (SOLID SURFACE) 
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• 200 

• 300 

• 400 

8500 

8600 

• 700 
8800 
8 900 
9000 
9100 
9200 
9300 
9400 
9500 
9600 
97 00 
9800 
9900 

10000 
10100 
10200 
10300 
10400 
10500 
10600 
10700 
10»00 
10900 
11000 
11100 
11200 
1 1300 
11400 
11500 
11600 
11700 
11900 
11900 
12000 
12100 
12200 
1 2300 
12400 
1 2500 
12600 
1 2700 
12800 
12900 
13000 
13100 
13200 
1 3300 
13400 
13500 
13600 
13700 
13800 
13900 
14000 
14100 
14200 
14300 
14400 
14500 
14600 
14700 
14800 
14900 
15000 
15100 
15200 
15300 
15400 
15500 
15600 
15700 
15800 
15900 
16000 
16100 
16200 


C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


12-APR-1981 23;29;54.94 
CHEMICAL REACTICN 4, CO +02 r C02 (SOLID SURFACE) 


ESI = ACTIVATION ENERGl FOR REACTICN 1 ( C A L/mOL E ) , E TC . 

Cl = PREEXPONEMIAL FACTOF (DIMENSIONAL) FOP REACTION 1, ETC. 

WHC, aCO = molecular NEJGHIS for CNhM, CCi 

OIS C heat of COMBUSTION FOR REACTION 1 (CAL/G CF CNHMi) 

02S = HEAT OF COMBUSTION FOR REACTICN 2 (CAL/G OF CC) 

OS s HEAT OF COMBUSTION FOR REACTION 3 (CAL/G OF CNH*-) 

XN = THE NUMBER OF CARBONS IN CNHM 
XM s the number of HYDROGEN IN CNHM 


OS = OIS ♦ 02S»XNArCC/WHC 

*+«++««» CATALYTIC BED PARAMETERS ♦♦♦♦♦»♦»»»»♦»» 

XLS = BED LENGTH (M) 

DS = CHANNEL HYDRAULIC DIAMETER (M) 

SOA = CIRCUMFERENTIAL LENGTH DIVIDED BY CHANNEL ARE A , = 4 /DS ( 1 /M ) 

AOAS = RATIO OF OPEN TO CLOSED AREAS 

CSS = HEAT CAPACITY OF SUBSTRATE (CAL/KG/K) 

ROS = DENSITY OF SUBSTRATE (KG/M3) 


READ (5,10) XPHl 

READ (5,10) IJ 

REAL (5,10) XFI 

READ (5,10) UIUS 

REAL (5,10) XLS 

READ (5,10) DS 

REAL (5,10) AGAE 

REAL (5,10) DX 

READ (5,10) DT 

READ (5,10) TMAX 

READ (5,8) IPRINT 

READ (5,8) N 

READ (5,8) NAS 

REAL (5,A) 01,02,03,04 

REAL (5,8) IT 

SOA = 4.0/DS 

CN = XN»wCC/wHC 

NPNAS = N + NAS 

OPERATING PARAMETERS 

TI = INITIAL, UPSTREAM (l.U.) TEMPERATURE (X) 
UI = l.U. VELCCITY (M/S) 

XPT s l.U. PRESSURE (ATM.) 

RGI = l.U. GAS DENSITY (KG/M3) 

YFI = l.U. MASS FRACTION OF HC 

YCI = l.U. MASS FRACTION OF CC 

YOI = l.U. MASS FRACTION OF 02 

XPhI S UPSTREAM HC/02 EQUIVALENCE RATIO 


UI = U1US»(1 .0 + , OASl/AOAS 

XMU = 0.000001458»SQRT(T1)/(1. 0+110. 4/TI) 

YFI = WHC4X^riI/29./((XN+XM/4.)44.76+XPHI) 
XLE*1 = 1.5» SORT (29./8HC) 

XLEN2 « 1.5» SORT (29./NCC) 

XLE123=XLEW144.667 

XLE223=XLEh2»».667 

XLOU = XLS/UI 

• + + * (;as property 4A************* 


CP = SPECIFIC HEAT OF GAS (CAL/KG/K) 

XKI a HEAT CONDUCTIVITY OF GAS AT l.U. (CAL/M/5EC/K ) 
ALPHA = l.U. GAS thermal DIFFUSIVITY (M2/SEC) 

GCaGAS CONSTANT FOR GAS MIXTURE (JOULE/KG/K) 

XLEhl a LEhlS NUMBER FOR CNHM IN MIXTURE 
XLEK2 a LEMIS NUMBER FOR CO IN MIXTURE 


XKI=.003+.000013*TI 
ROI a XPl/GC/TI 
ALPHA a XM/RCI/CP 
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16300 
16400 
16500 
16600 
16700 
16800 
16900 
17000 
17100 
17200 
17300 
17400 
17500 
17600 
17700 
17800 
17900 
18000 
18100 
18200 
1 8300 
18400 
18500 
11600 
11700 
11100 
1 8900 
19000 
19100 
19200 
19300 
19400 
19500 
19600 
19700 
19800 
19900 
20000 
20100 
20200 
20300 
20400 
20500 
20600 
20700 
20800 
20900 
21000 
21100 
21200 
21300 
21400 
21500 
21600 
21700 
21800 
21900 
22000 
22100 
22200 
22300 
22400 
22500 
22600 
22700 
22800 
22900 
23000 
23100 
23200 
23300 
23400 
23500 
23600 
23700 
23100 
23900 
24000 
24100 
24200 
24300 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 


c 

c 

c 

c 


c 

c 

c 


RFY = U1US*DS»HCI/XMU 


TO DESCRIBE CHE81CAL REACTION RATE EXPRESSIONS 
reaction 1 (GAS PHASE) 

D(CNH8)/DT tSS X P * » . 3 ♦ T» SORT ( V HC ) ♦ 1 0 2/EX P ( E 1 / R I) 
REACTION 2 (GAS PHASE) 

D(CC)/DT $tS YCC*SORT(YC2)*SCPT(yh2C)/EXP(E2/Rl ) 

REACTION 3 (SOLID SURfACE) 

D(CNHM)/DT S$S YHCS/EXP(E3/RT) 

REACTION 4 (SOLID SURFACE) 

D(CC)/DT $8$ YCCS/EXPCE4/RT) 

TO EVALUATE KINETIC RATE CONSTANTS 

BISS r C14S0RT (KHC)/(»C2/S0RT(1000. ) 

BIS = eiSS»XPI4».3*RCI*S0RT(RCI)8TI4SQRT(YFl)»Y01 
B1 = XL0U»B1S/RCI/YE1 
SSN02=S0RT(S0RT ( W02) ) 

B2S5 X C2/1000.»».75/5S802/SORT(WH20) 

SSRCI = SCiRT(SQRT(R0I)) 

SSYCIxS0RT(S0RT(Y0I)) 
B2S=B2SS8R0I*RCI/SSRCI»YE14SSY0I 
P2 = B2S8XL0U/RC1/YFI 

B3 X C38DS/ALPHA/NHC/3.66 
B4 X C48DS/ALPHA/WC0/3.66 

01 X GlS/CP/Tl 

02 X 02S/CP/TI 
0 X QS/CP/TI 
TAB X Yfl»0+l,0 
TABDM X TAB*TI 

El X ESl/II/1 .987 
E2 X ES2/T1/1.987 
E3 X ES3/TI/1.987 
E4 X ES4/II/1.967 
AD2UL X (ALPHA/CS/DS)/(UI/XLS) 


WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

■WRITE 

WRITE 


XPHI 


(6,21) N,NAS 
(6,20) TWAX 
(6,15) YCl, YEI, 

(6.71) I'US 

(6.72) UI 

(6,17) 01, 02, DX, DT 

(6.402) TI,XPR,R£Y, AD2UL, AOAS 

(6.409) TAB 

(6.410) TABDM 

(6.403) ES1, ES2, ES3, ES4 

(6.47) E1 , E2, E3, £4 
(6,35) Cl, C2, C3, C4 

(6.48) GC, XLEwl, XLEW2, XLS, OS 

(6.49) RCl, XKI, CP, ALPHA 
(6,401) XPl 

(6.43) XLCO 

(6.44) BIS, B2S 

(6.45) ei, B2, B3, B4 
(6,19) 

(6.24) 

(6.25) 

(6.24) 

(6,33) 

(6.24) (YC5(I), lxl,N,IT) 


(TS (I), Ixl,N,IT) 
(YFS(I), Ixl,N,IT) 


$$$$t$****$*$$*i1‘****$*$***$$ 

UPSTREAM CONDITICNS 

YF(1) X 1. 

YC(l) X 0.0 

yo(i) X 1. 

Td) X 1. 

R0(1) X 1. 

XP(1) X 1,0 
U(l) X 1. 

lt*$t$t*t*******$***$$$******$’t*¥t*$$$*****$* 


PAGE 3 
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24400 

24500 

24600 

24700 

24800 

24900 

24950 

25000 

25100 

25200 

25300 

25400 

25500 

25600 

25700 

25800 

25900 

26000 

26100 

26200 

26300 

26400 

26500 

26600 

26700 

26100 

26900 

27000 

27100 

27200 

27300 

27400 

27500 

27600 

27700 

27800 

27900 

28000 

28100 

28150 

28300 

28400 

28500 

28600 

28700 

28800 

28900 

29000 

29100 

29200 

29300 

29400 

29500 

29600 

29700 

29800 

29900 

30000 

30100 

30200 

30300 

30400 

30500 

30600 

30700 

30800 

30900 

31000 

31100 

31200 

31300 

31400 

31500 

31600 

31700 

31100 

31900 

32000 

32100 

32200 

32300 


START OF INTEGRATION LOOP 
KOUNT = 1 

XNUO = .8854SGRT(A0AS)4faCiRl(XPR4RFY) 

KOCH = 0.0 
R00T2 = 1.0 
555 CONIINUE 

ROOT = .5»(RC0T1 ♦ RCCT2) 

IFCFOOT.LE.O. nYNU=.0444/RD0T+3.46“l .344RCCT 

IF(ROOT.GT.0.1)YNU=0.011/R00T+3.66 

YKOCTsYNU-XNUO 

IF (ABSlYROCn.LT.0.01) GO TC 200 
IF (YKOCT.GT.0.0) ROOTlsRCOT 
IF (YROOT.LE.O.O) R00T2=P0CT 
GO TO 555 
200 CONTINUE 

DNASs=RC(i )*um 

INTEGRATION OF CUASl-STEADY GAS PHASE EQUATIONS FOR 
YF CFUEL MASS FRACTION) AND T (TEMPERATURE) IN DOWNSTREAM 
DIRECTION USING FOURTH-ORDER RUNGE-KUTTA SCHEME 
1 IS THE INDEX FOR DCwNSTREAM DISTANCE X 


C 

C 

c 

c 

c 

c 

c 

c 


DO 201 1=1, NFNAS 

1F(1.GT.N)DMAES=R0(1)4U(1)»ACAS/(ACAS+1.) 

4444444444444444444444444444444444444444W 

EVALUATE XJH AND XJD AS A FUNCTION CF X 

NONCIMENSIONAL HEAT CONDUCTIVITY 

*****«;4*44««4444444l4»»44**U4444«4*4«44 

IF(I.GT.l) GO TC 180 
IF(I.EO.l) X=0.0 
GO TO 181 

180 CONTINUE 
1X=I+1 
XI=IX-1 

X = X14DX 

181 CONTINUE 

XST = X»2.0»AD2UL»XK+ROOT 

IF (XST. LE. 0.1) XNU = 0.0444/'XET + 3.46-l .344XST 

IF (XST.GT.O.lj XNU=0.011/XST +3.66 

RNKCI)=XNU»XK/3.66 

XJH(I) = 4.04XNU 4AD2UL4XK 

IF CI.GT.N) XJHCI) = 0.0 

XJCKI) = XJh(l) 4XLE123 

XJC2CI) s XJH(I) 4XLE223 

W»«W**«»****»*«*4W«*«*W»««»4«»*«*4***44 

132 CONTINUE 
H=1 

DO 120 L=l,4 
IF (M.NE.l) GO ,0 102 
YFF = YF(I) 

YCC = YC( 1 ^ 

TT = TCI ) 

M = M + 1 
GO TO 119 

102 CONTINUE 

IF (M.NE.2) GO TO 103 
YFF=yF(I)+DX»XHY ( 1 )/2. 

YCC = YCCI) ♦DX*XKC(l)/2. 
TT=T(I)4XKT(l)4CX/2. 

M = M + 1 
GO TO 119 

103 CONTINUE 
IF(M.NE.3) GO TO 104 

YFF = YFCI) ♦ CX4XKY(2)/2. 

YCC = YCCI) + CX4XKCC2)/2. 

TT=T(I)+DX»XKT(2)/2. 

M = M+1 
GO TO 119 

104 CONTINUE 

YFF = YF(1) + CX»XKY(3) 

YCC = YC(I) + DX4XKC(3) 

TT=TC1) +DX4XKTC3) 
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IGNITION. F0h;17 


12-APR-1981 23:29:54.94 


PAGE 5 


32400 

32500 

32600 

32650 

32700 

32100 

32900 

33000 

33100 

33200 

33300 

33400 

33500 

33600 

33700 

33900 

33900 

34000 

34100 

34200 

34300 

34400 

34500 

34600 

34700 

34100 

34900 

35000 

35100 

35200 

35300 

35400 

35500 

35600 

35700 

35100 

35900 

36000 

36100 

36200 

36300 

36400 

36500 

36600 

36700 

36800 

36900 

37000 

37100 

37200 

37300 

37400 

37500 

37600 

37700 

37800 

37900 

38000 

38100 

38200 

38300 

38400 

38500 

38600 

38700 

38800 

38900 

39000 

39100 

39200 

39300 

39400 

39500 

39600 

39700 

39800 

39900 

40000 

40100 

40200 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 


119 CCMINUE 

' YH2C = MASS FRACTION OF H20 

TR2CCI) = .5»XM4WH2C»YF I» ( VF ( 1 ) - V F C I ) I / » HC 
IFdFF.LT.O.O) YFF=0.0 

! Ml = 61* XP ( I )»». 3»PC( I )4S0RT(R0(I ) )»T( 1 ) »SURT( YFF ) ♦ YC( I ) 
1/EXF(E1/TT) 

S6YC=SUR1 (SORT (YU( 1 ) n 

R017b=R0(I)*R0(l)/S0Rl(S0RT(KC(I))) 

t*2=fc2*R0175»YCC4SSYC»S0RTCYH2C(I))/EXP(E2/TT) 

XFYCU = -XJDim »(YFF- YF5 ( I ) ) /DM ASS -«1/DMAS5 
/XKC(L) = -XJD2CI) ♦CYCC-YCS(l))/ CMASS +Cn»m 1/DMASS -M2/DMA5S 
/ XKT(L) =>XJhCl)*(TT-TS(I))/DMASS+01»YFI*wl/DMASS+02»YF14W2/DMASS 

120 CONTINUE 

YFCl+l) =YF(I) +DX4(XKY (l)/6. +XKY (2)/3. +XKY (3)/3.+XKY (4)/6.) 
IF ( VF (I + l ) .LT.0.0 ) YF(1 + 1) =C.0 

YC(I + l) = YC(I) + CX*(XKC(n/6.-fXI<C(2J/3. + XKC(3)/3.4XKC(4)/6.) 

IF (YCCI + n.LT.O.O) YC(1 + 1) sQ.O 

T (I + l) =T (I) 4DX4CXKT m/6. +XKT (2)/3. +XKT (3)/3. + XKT (4)/6.) 
IF (T ll + l ) .LT .0.0) GC TO 601 
ADD=YC(I+1 )»YFI»M02/YCI/MC0/2. 

Y0(l+I)*1.0-XFH1»C1.0-YF(I+1))+ADC 

XP(I+1)=XP(I) 

R0(I+1)=XP(I+1) /T(I+1) 

UCI+l) sDHA5S/T(l+l ) 

EFF(I)=1 .-YF(I)-YC(I)402/0 
EFFI(I)s(T(I)-T(l))/(TAB-T(l)) 

201 CONTINUE 

END OF ONE RUNGE-KUTTA INTEGRATION CYCLE FOR QUASI-STEADY 
GAS PHASE EQUATIONS 


INTEGRATION OF UNSTEADY SOLID HEAT TRANSFER EOUATICF 
FOR TS (SOLID TEMPERATURE) IN TIME USING FCUPTH-QRDER 
RUNGE-KUTTA SCHEME 


C 

c 

c 

c 

c 

c 


DC 301 J =1 ,N 
MJsl 
DC 320 LJ=1,4 
IF (MJ.NE.l) GC TO 302 
TSS =TS(J) 

MJsRJ+1 
GO TO 319 

302 CONTINUE 
IF(MJ.Nt.2) GO TO 303 

TSS = TS(J) 4 DT*XKTS(l)/2. 

MJ=MJ+1 

GO TO 319 

303 CONTINUE 
IF(MJ.NE.3) go to 304 
TSS=TS( J)+DT4XKTS(2) /2. 

MJsMJ+1 

GO TO 319 

304 CONTINUE 

TSS = 1S(U) + DT» XKTS(3) 

319 CONTINUE 

XKTSCLJ)=RNK(J)»(T(J)-TSS) +B 3 * YF I ♦Q4R0 ( J ) » YFS ( J ) /EXP ( E3/TSS ) 

1 +f4*YFI»U24R0(J)»YCS(J)/EXP(E4/TSS) 

320 CONTINUE 

TS(J) =T6(J) ♦CT*(XKTS(l)/6.+XKT5(2)/3.+XKTS(3)/3.+XKTS(4)/b. ) 
YFS(J) = YF(J) /(l.+E34RO(J)/XLE123/RNK(J)/EXP(E3/TS(J))) 
YCS(J) = YC(J) /(l.+B4»R0( J)/XLE223/RNK ( J)/EXP(E4/TS( J) ) ) 

301 CONTINUE 

END OF ONE TIME INTEGRATION STEP 

LOOP FOR MARCHING FOPMARD IN TIME 

KOUNT = KOUNT+1 
KOU = KOUNT-1 
XKOUNT = KOUNT 
TIMEs(XKOUNT-l . )*DT 
IF (TIME.GT. TMAX) GG TO 999 
IF (MOD(KOU,IPRINT).NE.O) GC TO 200 
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IGNITION ,FOK ; 1 7 


40300 

40400 

40500 

40600 

40700 

40100 

40900 

41000 

41100 

41200 

41300 

41400 

41500 

41600 

41700 

41100 

41900 

42000 

42100 

42200 

42300 

42400 C 

42500 C 

42600 C 

42700 C 

42100 601 

42900 

43000 

43100 999 

43200 

43300 


12-APR-1981 23:29:54.94 
TIKE ,KCUNT 
(T (I),1=1,NENAS,1T) 

(TS (1), 1 = 1, N, IT) 

(YE ( 1 ) , 1=1 ,NPNAS, IT) 

( YC ( I ) , 1=1 ,NPNAS, IT) 
(YFSm,l = l ,N,U) 
(YCSC1),I=1,N,IT) 
(XJH(l),Isl,N,lT) 
(EFF(1),1=1,NPNAS,IT) 
(EFFT(I),1=1,NPNA5,1T) 


PAGE b 


WKITE (6,22) 

NRITE (6,13) 

WPnE(6,24) 

WRITE (6,9) 

WRITE(b,24) 

WRITE(6, 14) 

WRI1E(6,24) 

>»RITE(6,41 ) 

WRllE(b,24) 

WRITECb, 1 1 ) 

WRIIt(6,24) 

WRITE(6,42) 

WR1TE(6,24) 

WPITE(6,26) 

WR1IE(6,24) 

WRI1E(6,73) 

WRITE(6,24) 

WPIIE(6,74) 

WPIIE(6,24) 

WRnE(b,23) 

CO TO 200 


END OF WHOLE INTEGRATION LOOP 


CONTINUE 

WRITE (6,404) 

WR1TE(6,405) 

CONTINUE 

STCF 

END 
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BLOCK.. F0R;31 


12-APH-1981 23:32:53.58 


PAGE 1 


100 

200 

300 

400 

600 

700 

800 

900 

1100 

1200 

1400 

1500 

1600 

1700 

1100 

2000 

2100 


BLOCK DATA 

DIMENSION TS( 300),YFS( 300) ,YCS(300) 

COMMON /ACTI/ESl,tS2,ES3,ES4 
COMMON /MOLED/ ohC, WCC, m 02, MH20 
COMMON /THl/AN, XM, OIS, 02S 
common /TH2/ CP, GC, XPR 
COMMON /UPST/ YCI, YOl 
COMMON /1M/T5,YFS,YCS 

DATA ESI ,tS2,ES3 ,ES4/24 000. ,40000. , lOOoO. 
DATA aHC,InC0,*i02,aH 2C/44.,28.,32.,18./ 
DATA XN,XM,QlS,C2S/3., 8., 6491000., 24 15000 
DATA CP,GC,XPK/285. ,2.l312E-3,0.7/ 

DATA YCI, YOl/O. ,0,23/ 

DATA TS,yFS,YCS/3004].,30040.,30040./ 

END 


,17800./ 

./ 
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Example calculation 


* RUM IGNITIDN.exe; 10 
. IS 
1 0 0 0 . 

10 . 

. 1 

.0016 

Toi 
. 01 
1 . 0 0 1 
£ 0 
46 

so 

1.5E5 0.71E14 E.5E3 1 . E5 



MG' DF GRID PDIMT = 46 






r 4 

TI 

xpi 
UlUS 
X LS 
V3 

AOAS 
DX 
DT 

TMAX 
XPRIM7 

NAS 

Cl, C.2, C3 
IT 


C4 


80 


TMFiX = 1.00100005 

YDUJPSTjIMD = 0. £30000 YF aiPST? IMI> = 0. 0113888 EQUI RflTID = 0. 18000 

0007 


ftPPRGFiCH VEL = 10. 0000 

CHFIMMEL VEL a.U.> = 15.00000 
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01= ££. 77544 i!i£ Q£= 8.4736843 IiX= 0. 01 CnjOCiO IiT= 0. 

0 1 ij 0 0 (1 0 


TI» PRj REV? RDclUL* RORS = 1000. 00000 0. 70000 40S.31E38 0.13787 £ 

. 00000 

RDIRERTIC ELRME TEMP = 1. 44368c.' 

RDIRB FLRME TEMP <:K> = 1443.6£185 

E S !« S = 0 . £ 4 0 0 0 0 0 0 E + 0 5 0 . 4 0 0 0 0 0 0 0 E + 0 5 0. 10 0 0 0 0 0 0 E + 0 5 0.17800 0 

OOE+05 

EjS = 1£. 078510 £0. 130850 5. 03£71£ 6.85£'££8 

CiS = 0. 1500000000E + 06 . 0. 71 00000371E-t-14 0. £500000000E+04 0.1 00000 

OOOOE+06 


i3Cj XLEWl* XLEWc'j XLS> US = 0. 00£831 1. £17766 1.5E6551 0.100000 0 

.001600 

RDIj XKIk CP> RLPHR = 1.0586c'! 0. 016000 £85. 000000 0.5£98£E-04 


>iP I 


3. 00000 


L U = 0. 0066667 

FIS'. B£S. = 0.36601338E+05 0. 34534906E + 09 
Fl^ F£i. BSi. B4 = 0.£0£19805E+05 0. 19078£40E-*'09 0. 46881 531 E-*- 03 0. £946e391E-t-05 

INITIRL CDMIUTIONS CTS^ 

O.lOOOOE+01 O.lOOOOE+01 O.lOOOOE+01 O.lOOOOE+01 O.lOOOOE+01 O.IOOOOE+Ol O.IOOOOE 
+01 0. lOOOOE+01 0. lOOOOE+01 O.lOOOOE+01 

IMITIRL CDMDITIDMS (YPS,-. 

O.OOOOOE+00 O.OOOOOE+00 0. OOOOOE+00 0. OOOOOE+00 0. OOOOOE+00 0. OOOOOE+00 0. OOOOOE 
+00 0. OOOOOE+00 0. OOOOOE+00 O.OOOOOE+OO 

IMITIRL CQMDITIDNS: (YCS) 

O.OOOOOE+OO O.OOOOOE+OO O.OOOOOE+OO O.OOOOOE+OO O.OOOOOE+OO O.OOOOOE+OO 0. OOOOOE 
+00 O.OOOOOE+OO O.OOOOOE+OO O.OOOOOE+OO 
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T1ME = 


0 . c (1 0 0 U ij 0 0 


k:DUNT= 


gr: tempeprtlire a:> 

0. 1 OOOOE + Ol 0. 1 Oc.'?9E + CH 0.10351E+01 0. 1 CI393E+01 0. 1 04E'ijE + 01 0. 1 Ci437E+Ul U.1 0447E 
+01 0.10451E+01 0.10451E+01 0.10447E+01 

0. 1 04«jc'E+01 0. 10481E+01 0.10501E + 01 0.10521E+01 0.10541E+01 0. 1 05Ec'E+01 0. 1 0583E 
+ 0 1 0.1 08 04E+ 0 1 0.1 08E5E+ 01 0.1 0647E+ 0 1 

0. 1 0bfc.9E+01 0.10891E+01 0.10714E+01 0.10737E+01 0.10760E+01 0.10784E+01 
SDL in TEMPERRTLIRE <;TS> 

0.11390E+01 0.10788E+01 0.10631E+01 0. 10548E+01 0.10477E+01 0.10487E+01 0. 1 0384E 
+ 01 0. 10349E + 01 0. 1 03c.'lE+01 0.10295E+01 

CNHM MOSS FPfiCTIDH (YF> 

O.lOOOOE+01 0.8£811E+00 0.7£848E+00 0.64850E+00 0.58054E+00 0.5£141E+00 0.4893£E 
+00 0.4£307E+00 0.38150E+00 0.34400E+00 

0.33115E+00 0.3£383E+00 0.31648E+00 0.3090£E+00 0.30153E+00 0.£9398E+00 0.£8838E 
+ 00 0.£7871E+0o 0.£7100E + U0 0. lci‘83*£4E+00 

0.£554£E+00 0.£4756E+00 0.£3966E+00 0.£3171E+00 0.££373E+00 0.£1571E+00 
YCD = 

O.OOOOOE+00 0.1 1484E-01 0.££138E-01 0. 31487E-01 0. 395£c'E-01 0.46303E-01 0.519£3E 
-01 0.58483E-01 0.60031E-01 0.6£681E-01 

0.73834E-01 0.87881E-01 0.10073E+00 0.11418E+00 0. l£76£E+00 0.14103E+00 0. 15440E 
+00 0.1877£E+00 0. 18098E+00 0. 19417E+00 

0.£07£6E+00 0.££0£4E+00 0.£3309E+00 0.£4579E+00 0.£5833E+00 0.£7088E+00 
SUPFRCE HC MftJS FRRCTIDM CYFS) 

0.53851E+00 0.£4678E+00 0.£0088E+00 0.173£6E+00 0. 15£71E+00 0. 13598E+00 0.1E175E 
+00 0.10985E+00 0.99795E-01 0.90590E-01 


YCO<:SLIRFF(C:E> = 


O.OOOOOE+00 0.£84£9E-0£ 0.48789E-0£ 0.88753E-0£ 0.8641£E-0£ 0.10176E-01 0.11484E 
-01 0.1£6£8E-01 0.13836E-0'1 0. 14440E-01 


XJH 

0.11877E+0£ 0.3££71E+01 0.£6856E+01 0.£4'.99E+01 0.£33£8E + 01 0.££489E + 01 0.£1875E 
+01 0.£1536E+01 0.£1484E+01 0.£1403E+01 

EFF I C I EMC Y (CFtRBDM-BRLRMCE) 

O.OOOOOE+00 0.16939E+00 0.£6670E+00 0.34465E+00 0.41088E+00 0.4885£E+00 0.51939E 
+00 0.56464E+00 0.60544E+00 0.64£37E+00 

0.85£79E+00 0.65718E+00 0.66163E+00 0.66614E+00 0.67071E+00 0.87534E+00 0.88004E 
+00 0.68480E+00 0.68963E+00 0.69458E+00 

0.69949E+00 0.70453E+00 0.70984E+00 0.7148£E+00 0.7£008E+00 0.7£541E+00 
EFFICIENCY CTEMPERRTURE? 

0. OOOOOE+00 0.6£838E-01 0.79086E-01 0.88607E-01 0.94869E-01 0.98477E-01 0. 1 0086E 
+00 0.10160E+00 0.10161E+00 0.10085E+00 

0.10409E+00 0.10848E+00 0.11£9£E+00 0.11743E+00 0. l££00E+00 0. 1£663E+00 0. 13133E 
+00 0.13809E+00 0.1409£E+00 0.14588E+00 

0.15079E+00 0.1558£E+00 0. 16093E+00 0. 16611E+00 0.17137E+00 0.17870E+00 
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TIME = 


0 . 4 Cl Cl Cl u Cl Cl 1 


KDUMT= 


41 


i3l=i:i TEMFERftTIJRE CT;' 

Cl. 1 OOiKiE+dl 0. 10491E+CH 0. 1 CI635E+CU Ci. 1 Cj7£fjE+rjl 0. 1 CI776E+CH O. 1 Ci81 lE + CU Ci.lCi833E 
+ 01 0. 1 0843E + CU Cl. 1 Oe45E + Cil 0. 1 Ci841E + Cil 

0. 1 0861E + 01 Cl. 1CI8S9E + CI1 0. 1 0917E + CU 0. 1 0945E+CU 0. 1 CI975E + CU O. 1 1 0Ci4E + i:il O.llCiSSE 

+ Cil Cl. 1 1 085E+CI1 0. 1 1 CI97E + 01 0. 1 1 lc'9E+CU 

Cl. llie.lE+01 0. 1 i 194E + CU O. 1 l££8E+i;U Ci.ll£6£E+01 0. 1 1£97E+CU 0.1133£E+CH 

SOLID TEMPERATURE CTS:> 

0. 1££13E+CU 0. 1 1518E + CU 0. 1 1£46E+Cil 0. 1 1 CI83E+CU 0. 1 0964E + CU 0. 1 Ci869E+ CU Ci.lCi789E 

+ 01 0. 1 Ci7££E + 01 0. 1 0667E+01 0.10817E+01 

CNHM MASS FRACTIDM CYF> 

0. 1 OCiOOE + Cil 0. 80756E+0CI 0. 701 37E+00 0. fel707E+Ci0 0. 548CilE+00 0. 48487E + CI0 CI.43108E 
+ 00 0. 38389E + 0CI 0. 34185E+0CI CI.30433E+00 

0. £89£5E+0Ci 0. £7908E+0Ci 0. £6883E+00 0.£5851E+00 0. £4813E+CiCi 0. £3789E+CiCi 0.££719E 
+ 00 0.£1866E + 00 0. £0809E+CiCi 0.195508 + 00 

CI.16489E + 00 0.174£9E+00 0.183708+00 0. 15314E+00 0. 14£84E+00 0. 1 3££1E+ 0 0 
YC:0 = 

0. OOCiOOE + OCi 0.13£££E-01 0. £877£E-01 0.39£34E-01 0.50££OE-01 0.598088-01 0.87399E 

-01 0.73857E-01 0. 78404E-01 0.81781E-01 

0. 98870E-CU 0. 1 1495E+00 0. 13388E + 00 0. 150808+00 0.186098+00 0. 185£9E + CiO 0.£0£18E 
+ 00 0. £1 889E + 00 0. £3479E + 00 0. £5Ci4CiE + 0Ci 

0. £85488 + 0 0 0. £7995E + 00 0. £93788+00 0. 3Ci86£E+0Ci 0.319078+00 0.33 0438+00 
SURFACE HC MASS FRACTIDM CYFS) 

0. 484£4E+00 0.198388+00 0.184748+00 0. 143808 + 00 0. l£76£E + 0Ci 0. 1 14 038+00 0. 1 0318E 
+00 0.9££77E-01 0.635108-01 0.75484E-01 

ycocsurface;. = 


Cl. Cl Cl Cl Cl CiE+ Cl Cl 0.£Ci35£E — 0£ 0. 4£Ci4£E— Ci£ 0. 8‘3i778E— 0£ 0.845718— 0£ 0. 1031788—01 0. 180908 
-01 0.138818-01 0.150848-01 0.181408-01 

XJH 

0.118778+08 0.38884F+01 0.87E97E+01 0.850958+01 0.838588+01 0. £30398+01 0.884358 
+01 0. ££1808+01 0. ££1 1£E+01 0. ££0498+01 

EFF I C I EMC Y CCARUDN-BALAMCE) 

0.000008+00 0.189548+00 0. £98818+00 0.374408+00 0.443088+00 0. 50838E+CI0 0.554858 
+00 0.800088+00 0.841098+00 0.877888+00 

0. 889688+00 0.69598E+00 0.708878+00 0.708738+00 0.715318+00 0.788008+00 0.788888 
+00 0.735778+00 0.748848+00 0.750038+00 

0.757368+00 0.764818+00 0.778408+00 0.780118+00 0.787958+00 0.795918+00 
EFFICIENCY (TEMPERATURE) 

0.000008+00 0.110748+00 0.143048+00 0.188338+00 0.174858+00 0.188988+00 0.187758 
+00 0.190108+00 0.190588+00 0.189538+00 

0.194078+00 0.800318+00 0.806668+00 0.813188+00 0.819708+00 0.886408+00 0.833888 
+00 0.840168+00 0.847838+00 0.854438+00 

0.881758+00 0, £69818+00 0.876798+00 0.884518+00 0.898348+00 0.300318+00 
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TIME = 


0 • 59999996 


K,DUNT= 
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GftS TEMPEPfiTLIRE a> 

O.lOOOOE+01 0. 1 0637E+01 0. 1 0848E+01 0. 10979E+01 0. 1 1 067E-t-01 0.111E7E + 01 iXlllbbE 
+ i;H CI.11188E + 01 0. 11198E+01 O.U197E+i:n 

O. 1 1££5E+01 0. 1 1£6£E+01 0. 1 1301E+01 0. 1 r34CiE+01 0. llC-sSOE + CH 0. 1 14£1E+01 0. 11463E 
+ 01 0.115 06E+ 0 1 0.11 55 0E+ 0 1 0 . 1 1 594E+ 0 1 

0.11840E+01 0.11686E+01 0.11733E+01 0.11780E+01 0. U8£8E+01 0.11876E+01 
SOLID TEMPERfiTUPE CTS> 

0.1£635E+01 0. 1£134E+01 0.11791E+01 0.11578E+01 0.11419E+01 0.11£89E+01 0.11179E 
+01 0.11088E+01 O.llOllE+01 0.10940E+01 

CHHM MRSS FRPlCTION (YF> 

O.lOOOOE+01 0.79689E+00 0.68495E+00 0.59669E+00 0.5££39E+00 0.45843E+00 0.40£78E 
+00 0.35401E+00 0.31086E+00 0.£7£70E+00 

0.£5498E+00 0.£4161E+00 0.££817E+00 0.81489E+00 0.£0118E+00 0. 18787E+00 0. 17418E 
+00 0.16074E+00 0. 14739E+00 0. 13416E+00 

0.1£109E+00 0.108£5E+00 0.95689E-01 0.8347£E-01 0.71677E-01 0.60388E-01 
YCO = 

0. OOOOOE+00 0. 14469E-01 0.30578E-01 0.46068E-01 0.60086E-01 0.7££38E-01 0.8£367E 
-01 0.90438E-01 0.96435E-01 0.10054E+00 

0.11961E+00 0.14££1E+00 0. 164££E+00 0. 18553E+00 0.£060£E+00 0.££560E+00 0.£4413E 
+ LIU 0.£E'148E+O0 O.£775OE+0U 0.£9c'06E+00 

0.30499E+00 0.31614E+00 0.3£534E+00 0.33£44E+00 0.337£8E+00 0.3397£E+00 
SURFACE HC HASS FRACTION <:YFS> 

0.430££E+00 0. 1675£E+00 0. 14035E+00 0. l£3£3E+00 0.10988E+00 0.97974E-01 0.87588E 
-01 0.78891E-01 0.709E1E-01 0.63536E-01 

YCO (SURFACE) = 

O.OOuOOE+00 0. 181 0£E-0£ 0.35915E-0£ 0.57604E-0E' 0.79705E-0E' O.lOlllE-01 0. 1£099E 
-01 0.13980E-01 0.15616E-01 0. 16983E-01 


;cjH 

0.11877E+08 0.38865E+01 0.£76£8E+01 0.£5486E+01 0.£4£91E+01 0.£3501E+01 0.££915E 
+ 01 0.££748E + 01 0. c'c'ES^tjE + Ol 0.££E-33E + 01 

EFFICIENCY (CARBON-BALANCE) 

O.OOOOOE+00 0.£0056E+00 0.30840E+00 0.393£9E+00 0.46453E+00 0.58586E+00 0.57930E 
+ 00 0.6£6S-!£E+00 0.66817E+00 0.70543E+00 

0.71900E+00 0.7£746E+00 0.73610E+00 0.74495E+00 0.75400E+00 0.763E5E+00 0.77E71E 
+00 0.78238E+00 0.79E25E+00 0.80E31E+00 

0.81E56E+00 0.8££98E+00 0.83354E+00 0.84421E+00 0.85495E+00 0.86.571E+00 
EFFICIENCY (TEMPERATURE) 

O.OOOOOE+00 0.14356E+00 0.19111E+00 0.££064E+00 0.84056E+00 0.£5404E+00 0.£6£75E 
+00 0.26777E+00 0.26994E+00 0.26975E+00 

0.27609E+00 0.28455E+00 0.29320E+00 0.30204E+00 0.31109E+00 0.32034E+00 0.32980E 
+00 0.33947E+00 0.34934E+00 0.35940E+00 

0.36965E+00 0.38007E+00 0.39063E+00 0.40130E+00 0.41204E+00 0.42280E+00 
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TIME= 


Cl . 8 0 Cl Cl Cl Cl Cl 1 


KOUHT= 


81 


i3RS TEMPERfHTUPE CT;' 

0. 1 OCiOCiE+CU Cl. 1 Ci?38E+CH 0. 1 1 OC 18 E+OI 0.11ie3E+CU Ci. 11 3Cit.E+Cil Ci. 1 139£E + Cil Ci.ll45£E 
+ 01 i;M1490E+i;U CMISIEE+OI 0,11519E+01 

0. 1 lfi57E+01 0. I IfeOEE+Ol 0. 1 1&56E+U1 0. I l?i:i8E + CU 0. 1 1781E+CU 0. 1 1816E+CU Ci.ll871E 
+01 0.119£8E+01 0.11988E+01 0. 18045E+01 

Cl. 1 c* 1 04E+ 0 1 0.1£lb£E + i'il 0. l£££CiE+01 C*. 1££76E+01 0. 1£'330E+01 0. 1£380E+Cil 
SOLID TEMPERftTURE <JS> 

0. 1£84£E+01 0.1£637E+01 0. 1££63E+01 0. 1£0£1E+01 0.118368+01 0.11683E+01 O.ll^i^lE 
+ 01 0. 1144£E+Cil Cl. 1 1347E + CU 0.11£5eE+01 

CMHM MfiSS FRFiCTIOM <;YE> 

O.lOOOOE+01 0.78985E+00 0. 67457E+00 Ci.58£e6E+00 0.505478 + 00 0.438788+00 0. 3808EE 
+00 0.330138+00 0.885518+00 0.846368+00 

0.885748+00 0. 80896E + CiCi 0. 1988 08+ 00 0. 175498+00 0. 158888 + 00 0. 14£44E+CiO 0.1£6£5E 
+00 0.110388+00 0.949508-01 0.800698-01 

0.658798-01 0.585398-01 0.408888-01 0.891468-01 0.195118-01 0.115538-01 


VCD = 

0. OOOOOE+OO 0.153848-01 0.335748-01 0.518488-01 0.688088-01 0.837868-01 0.968188 
-01 0.106098+00 0.113888+00 0.117988+00 

0. 14070E + i:ii;i 0. 1671 18+00 0. 198048+00 0.815898+00 0.836598+00 0. £55718+00 0. £78388 
+00 0.886388+00 0.897888+00 0.304988+00 

0.309188+00 0.309708+00 0.306388+00 0.899138+00 0.887998+00 0.873068+00 
SURFftCE HC MftSS FR8CTI0M CYRS;' 

0. 414538+00 0. 14783E+CICI 0.183388+00 0. 108458+00 0.964888-01 0.859048-01 0.764188 
-01 0.684488-01 0.609598-01 0.540058-01 

YCOCSURFfiCE) = 

0.000008+00 0.13£53E-0£ 0.309878-08 0.516958-08 0.738368-08 0.960908-08 0.117378 
-01 0.137998-01 0.155998-01 0.171198-01 


YJH 

0. 1 16778+08 0. 33 03CIF KH 0.878768+01 0.857958 + 01 0.846468+01 0.838898+01 0.833888 
+ 01 0. £38438+01 CI.83806E + 01 0.83168E+01 

EFFICIEMCYCCFIRIiDM-FflLHNCE) 

0.000008+00 0.806888+00 0.318138+00 0.405868+00 0.479568+00 0.543018+00 0.598858 
+00 0.646798+00 0.689658+00 0.787988+00 

0.743668+00 0.754698+00 0.766038+00 0.777688+00 0.789658+00 0.801938+00 0.814508 
+00 0.887338+00 0.840388+00 0.853598+00 

0.866868+00 0.880098+00 0.893188+00 0.905788+00 0.917848+00 0.989058+00 
EFF I C I EMCY <TEnPERFlTURE> 

0.000008+00 0.166008+00 0.887818+00 0.866708+00 0.894348+00 0.313898+00 0.387318 
+00 0.335938+00 0.340788+00 0.348488+00 

0.351008+00 0.368038+00 0.373378+00 0.385038+00 0.397008+00 0.409888+00 0.481858 
+00 0.434688+00 0.447738+00 0.460938+00 

0.474818+00 0.487448+00 0.500478+00 0.513138+00 0.585198+00 0.536398+00 
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TIME= 


1 . 0 0 (1 0 0 0 0 Cl 


k:dunt= 


ICil 


GRS TEMF'ERRTLIRE 'CT'' 

Cl. 1 aCiOOE+01 0.10805E + CI1 d. 1 1 1E9E+CU 0. 1 1345E+CU 0. 1 l^CUE + CU 0. 1 lElEE + CU d. 11698E 

+ Cil d.ll755E + dl Cl. 1 179c'E-»-dl d.ueilE + Cil 

d.ll8ElE+Cil i;i.ll9£SE + dl d.ll986E + dl Cl. l£d5£E+dl Ci. l£l£dE + dl Ci.l£l&9E + dl Ci.l££59E 

+dl d.l£33dE+dl d.l£4dlE+dl d.i£47dE+dl 

Cl. 1E536E+CU d.l£597E+dl Ci. 1£E51E + CH Ci. 1EE.9E.E + CH 0. l£73CiE + Cil Ci.l£757E+dl 

SOL in TEMPERRTURE (TS:' 

Ci.l£94£E + dl d.l3Ci41E+dl Ci. 1£666E+CH Ci. 1£41£E+CU Ci. l££13E+dl Ci. l£Ci45E+dl CI.11899E 

+ Cil d.ll778E + Cil Cl. 1 167CiE + Cil d.ll568E + Cil 

CNHM MftSS FRACTION CYF> 

Cl. 1 ddCiOE+Cil d.7860£E-t-Cid Ci. 66774E-t-Cid 0.57313E+0CI 0. 49£91E+0Ci 0. 4£358E+CiCi Ci. 363£5E 
+ CiCi Cl. 31 CiSdE + Cid Cl. £64££E-*-CiCi Ci.££388E+CiCi 

Ci.£CiCi££E+Cid Cl. I8C1C1IE + C1C1 Cl. 15997E+CICI d.l4Ci£lE+CiCi Ci. l£Ci85E+CiCi Ci. 1 Ci£Ci5E+CiCi CI.83997E 
-dl Cl. 869CilE-dl d.51Ci£5E-dl Ci.36667E-dl 

Cl. £4183E-Cil CI.13877E-CI1 Ci.6187dE-0£ Ci. 14578E-Ci£ Ci. OdddCiE+CiCi Ci. CiCiCiCiCiE-»-CiCi 


YCD = 

Cl. dCiCiCiCiE + CiCi Cl. 159CI4E-CI1 CI.35886E-CH d.56611E-dl Ci.78315E-dl Ci. 93878E-CH Ci.lCi864E 
+ C1C1 Cl. l£Ci£CiE+CiCi Cl. l£84CiE-*-dCi d. 13347E + CiCi 

Cl . 1 59 1 EE-t- Cl Cl Cl . 1 88 Ci4E-*- Ci Ci Ci . £ 1 4 Ci3E+ Ci Ci 0 . E3868E+ Ci Ci 0 . £5558E+ 0 Ci 0 . £7 d£5E+ Ci Ci d . £8 Ci33E 
+ C1C1 Cl. £8548E-*-CiCi Ci. £8534E+CiCi Ci.£7984E+CiCi 

Cl. £6894E+CiCi Ci.£5£85E+CiCi Ci.£3195E + CiCi Ci. £d8e4E-*-CiCi Ci. 17831E+CiCi Ci. 15Ci73E+CiCi 
SURFACE HC MASS FRACTION CYFS) 

Cl. 4Ci7££E+CiCi Ci. 13348E + CiCi d.lll£7E+CiCi Ci.97579E-Cil Ci.8e.488E-Cil Ci.76871E-Cil Ci.6?731E 
-dl Ci.6Cil91E-Cil Ci.5£984E-dl Ci. 46£e5E-Cil 

YCDCSURFACE) = 

Cl. CiCiCiCiCiE-*-CiCi d.n34CiE-Ci£ Ci.£7£ldE-Ci£ Ci. 468Ci6E-Ci£ Ci.e.8d86E-0£ Ci.9d£87E-Ci£ Ci.lllAlE 
-dl Cl. 133£1E-CU d. 1518CIE-CI1 0. 16751E-CU 

;::jh 

Cl. 1 1677E+0£ Cl. 33 1 43E+ Cl t Ci. £8Cife3E'*’Cil d.£6d39E + dl Ci.£4936E+dl 0. £4£15E+dl 0.£3863E 
+ dl Cl. 83878E + CI1 Ci. £3C ..■.7E-*-CU O.E3841E+CU 

EFFICIENCY CCARBON-BALANCE) 

d.OdOOOE+OCi d.£lC5£E+CiCi 0.3£445E+0Ci d. 41456E-t-dCi Ci. 49Ci49E+Cid d.55600E+CiCi Ci,6131lE 
-t-OCi 0.66335E+0d 0.7Ci785E+dCi 0.7471dE+dCi 

0.76517E + 0CI 0.779CI9E+0CI d.79347E+0d d.8d830E+00 0. 8£355E+CiCi 0.83916E+0d Ci.8550£E 
■t-OCi 0.87100E+00 0.88690E-*-00 0,90£46E+00 

CI.91733E+0CI 0.9311EE+00 0.94335E+00 0. 95355E+0CI 0.96ieiE+0Ci Ci. 96721 E+ 00 
EFF 1 C I EHCY <TEMPERATURE> 

0. OOOOOE+00 0.18145E+00 0.25441E+00 0.30314E+00 0.33837E+00 0.36422E+00 0. 38283E 
+00 0.39570E+00 0.40395E+00 0.40818E+00 

0.41940E+CICI 0.4333EE+00 0.44770E+00 0.46253E+00 0.47778E+00 0.49339E+00 0.509E5E 
+00 0.52523E+00 0.541 13E+00 0.55669E+00 

0.57156E+00 0.58535E+00 0.59758E+00 0.60778E+00 0.61546E+00 0.62146E+00 
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APPENDIX II 


LISTING OF COMPUTER PROGRAMS FOR CHAPTER II 

The computer programs for Part II (steady state) are quite 
similar, in structure, to those in Part I. In addition to using a 
three-reactions scheme as discussed in the text of Part II, only the 
final results are printed out (no intermediate transient) . 

More Detailed Instruction (for CALYl. FOR; 75): 

1. Read-in data (lines 10400 - 11475): 

XPHI = (() = fuel-air equivalence ratio 

TI = T(0,0) = upstream, initial temperature (K) 

XPI =^(0,0) = upstream, initial pressure (Atm.) 

UIUS = ^'ref ~ 'upstream (outside of the bed) initial gas 
velocity (m/s) 

XLS = reference length (m) , suggested to use 0.1m (10 cm) 

DS = d* = channel hydraulic diameter (m) 

AOAS = A*/A^* = open-to-close area ratio 
DXX = nondimensional step size in X (space) inside the 
catalytic bed. It is nondimensionalized by XLS. 

So if XLS = 0.1 m and DX = 0.01, the dimensional 
step size is 0.001 m or 1mm. 

XAS = ratio of step sizes between downstream and inside-bed 

regions. For example, (XAS) x (DXX) is the integration 
step size in X - direction in the downstream space. 

DT = nondimensional time step size, recommended to be 0.2 0*1 

for steady calculations. 

TMAX = maximum nondimensional time allowed for computation. 

It should be large enough to allow the steady state 
to be reached. For bed length less than 5 cm, TMAX = 
8.001 is suggested, for bed length between 5-12 cm, 
TMAX = 12.001 is suggested. 



N = number of grid points (in X direction) in the 

catalytic bed. The nondimensional catalytic bed 
length is given by (N-1) x (hX) . For example, 
if N = 51, DX = .01, then the nondimensional bed 
length is equal to 0.5. If XLS = .Im, then the 
dimensional bed length is .05m (5 cm). 

NAS = nimber of grid points (in X - direction) in the 
downstream after-bed space. The length of the 
downstream distance can be calculated in a similar 
manner as in the above. 

IT = interval in X for printout. For example, if DXX =.01, 

IT= 5, the printout will be at X = 0, 0.05, 0.1, 0.15, 

etc .... 

CO, Cl, C2, C3, C4 = C^, C^, C^, C^ (See Table I) 

IPUTIN = 1 (always) 

Output data (lines 37000 -38586): 

The following quantities will be printed out as a 
function of downstream distance: 

Gas and solid temperatures, Emission Indices (xlOOO) for the 
original hydrocarbon, pyrolyzed hydrocarbon, total hydro- 
carbon and CO in the gas phase, the original hydrocarbon and 
CO at the wall and the combustion efficiency. 

If TMAX is large enough, the above printouts are the 
steady-state profiles. 
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FACIF 1 


100 
200 
300 
400 
500 
525 
537 
550 
575 
600 
• 00 
900 
1000 
1100 
1200 
1400 
1450 
1475 
1500 
1600 
1700 
1100 
1150 
1175 
1M7 
1900 
2000 
2100 
2200 
2300 
2400 
2500 
2600 
2700 
2*00 
2900 
3000 
3100 
3200 
3300 
3400 
3500 
3600 
3700 
3100 
3900 
4000 
4100 
4200 
4300 
4400 
4500 
4600 
4700 
4100 
4900 
5000 
5100 
5200 
5300 
5400 
5500 
5600 
5700 
5(00 
5150 
5900 
6000 
6100 
6200 
6300 
6400 
6500 
6600 
6650 
6675 
6617 
6743 
6799 
6100 
6900 


****t*^t***4*********4***tf******t**^t 

THREE GAS-PHASE REACTICN STEFS 

VARIABLE X STEF SIZE C EEC VS. AFTER-BED SPACE ) 

PARARFIER NX=300 

CQHHCN /ACTI/ ESO.tSl ,ES2,ES3,ES4 
COHFCN /MOLEG/ »hC, RCC, WC2, AH2G 
CORRCN /THl/XN, XM. QIS, 02S 
COMMON /TH2/ CP. GC, XPR 
COMMON /UP5T/ TCI, VOI 
COMMON /INl/TS.YFS.YCS.VES 
COMMON /B12345/P7A1,61A2,BTA3,ETA4,BTA5 
COMMON /GAMMA/GAMA 

dimension YF(NX), T(NX), TS(NX), YO(NX), YFS(NX) .RCCNX), U(NX) 
DIMENSION XKYC4), XKTC4), XKl S ( 4 ) . XKC ( 4 ) , YC(NX), YCS(NX) 
DIMENSION XJH(NX), XvlDKNX), XJ02(NX) , YH20(NX) .XPCNX) 
DIMENSION EFF(NX), RNK(NX) 

DIMENSION A(NX,3), XKE ( 4 ) , YE ( NX ) , YTHC ( NX ) , YES ( N X ) 

EQUIVALENCE ( A ( 1 , 1 J , TS ( 1 ) ) , ( A C 1 , 2 ) , YFS (1 5 ) , ( A ( 1 , 3 ) , YCS ( 1 ) ) 

7 FORMAT (4E20.10) 

• FORMAT (14) 

9 FORMAT (IHO.'SCLID TEMPERATURE (TS)'/) 

10 FORMAT (6F12.6) 

11 FORMAT (IHO, 'SURF ACE HC MASS FRACTION (YFS)'/) 

12 FORMAT (//) 

13 FORMAT (IHO, 'GAS TEMPERATURE (T)'/) 

14 FORMAT ( IHO, ’ORIGINAL CNhM (YF)'/) 

15 FOHMATdHO, 'YC(UFST,1NI) =', F10.6,bX, 'VF (UPST.INI) =',F10.7, 

1 bX.'tOUI RATIO s', F15.9//) 

17 FORMAT (lH0,'Ql=',F15.7,5X,'O2s', F 1 5 . 7 , 5X , ' DXs ' , F 15 . 7 , 5X , ' DTs ' , 
1F15.7//) 

II FORMAT (IHO, 'Ns' ,110//) 

19 FORMAT (IhO, 'INITIAL CCNDITICNS (TS)',/) 

20 FORMAT (1H0,'TMAX s ',M5.P) 

21 FORMAT (1H0,'NC. OF GRID POINT s', 2110//) 

22 format ( IHO, 'TIMES ' ,F15.M , 5X, ' MOUNTS ', II 0/ ) 

23 FORMAT (///) 

24 FORMAT (10E12.5) 

25 FORMAT (IhO, 'INITIAL CCNDITICNS (YFS)',/) 

26 FORMAT (IHO, 'XJH '/) 

31 FORMAT (IriO,' CC MASS FRACTION (YC)'/) 

32 FORMAT (IHO, 'SURFACE CO MASS FRACTION (YCS)'/) 

33 FORMAT (IHO, 'INITIAL CCNDITICNS (YCS)',/) 

- . > ES,S s • , 6F12.6/) 

' C,S s ' , 5E20.10/) 

' A,S s ' , 6F12.6/) 

' BE,S s ' , 6F12.6/) 

YCC s' ,/) 

YCC(SURFACE) *',/) 

L/U s', F15.7/) 

44 FORMAT (IHO,' BIS, B2S, s', 4E15.I/) 

45 FORMATdHO. ‘ bo-4 s'.SeIs.I/) 

" , 8 ) 

' E,S s',6F12.6/) 

'GC, XLENl, XLER2, XLS, DS »',6F12.6 /) 

'ROI, XKl, CP, ALPHA s', 3F1 2 . 6 , E 1 2 . 5/ ) 

'APPROACH VEL s «,F10.4/) 

'CHANNEL VEL (I.U.) s ',F10.5/) 


34 FORMAT (IHO, 

35 FORMAT (IHO, 
3b FORMAT (IHO, 
37 FORMAT (IHO, 

41 FORMAT (IHO, 

42 FORMAT (IHO, 

43 FORMAT (IHO, 


M3 rUKMWiljnU,' 

46 FORMAT (4E20 

47 FORMAT (IHO, 
41 FORMAT (IHO, 
49 FORMATdHO,' 

71 FORMAT (IHO, 

72 FORMAT (IHO 


401 FORMAT (IHO, 

402 FORMAT (IHO, 

403 FORMAT (IHO, 

404 FORMAT (IHO, 

405 FORMAT (IHO, 

406 FORMAT (IHO, 

407 FORMAT (IHO, 
401 FORMAT (IHO 


73 FORMATdHO, • EFFICIENCY (CARBON-BALANCE) ' ,/) 


XPI s' F20.S/) 

Tl, PR, REY, AD2UL, AOAS =',5F12.5/) 
EE,Ss’5E20.e/) 

GAS TEMPERATURE IS NEGATIVE') 

CNHM CONCENTRATICN IS NEGATIVE') 

CC CONCENTRATICN IS NEGATIVE') 

MU* runi-ff* vxnu,' YFF IS NEGATIVE') 

409 FORMATdHO, 'ADIABATIC FLAME TEMP s',F10,6/) 

410 FORMATdHO, 'ADIAE FLAME TEMP(K) =',F15.5/) 

411 FORMATdHO, 'BETA 1-5 s',5F10.4/) 

412 FORMATdHO, 'PYRCLYZED CNHM (YE)',/) 

413 FORMATdHO, 'TOTAL CNHM (YTHC)',/) 
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PAGE 2 


3 


7000 

7100 

7200 

7300 

7400 

7450 

7475 

7500 

7600 

7700 

7800 

7900 

8000 

1100 

1200 

1300 

• 400 

• 500 
8600 

• 700 
8800 

• 900 
9000 
9100 
9200 
9300 
9400 
9500 
9600 
9700 
9900 
9900 

10000 
10100 
10200 
10300 
10400 
10500 
10600 
10700 
10900 
10900 
11000 
11100 
11200 
11300 
11400 
11425 
11437 
11450 
11475 
11500 
11600 
11650 
11700 
11900 
11900 
12000 
12100 
12200 
12300 
12400 
12500 
12600 
12700 
12800 
12900 
1 3000 
13100 
13200 
13300 
13350 
13375 
13400 
13500 
13600 
13700 
13900 
13900 
14000 
14100 


***** INPUT PARAMETERS ***** 

***** CHEMICAL PARAMETERS »»«♦* 


CHEMICAL REACTION 0, 

CHEMICAL REACTION 1, 
CHEMICAL REACTION 2, 
CHEMICAL REACTICK 3, 
CHEMICAL REACTION 4, 


CNHM = CNHM (GAS PHASE PYPCLYEIS) 
(ISOTHERMAL REACTION) 

CNHM 4 C2 = CO +H2C (GAS PHASE ) 

CO + 02 = C02 (GAS PHASE) 

CNHM +C2 = C02 + H20 ( SOLIL SURFACE) 
CO + 02 = C02 (SOLID SURFACE) 


ESI s ACTIVATION ENERGY FOP REACTION 1 (CAL/MOLE), ETC. 

Cl = PREEXPONENTIAL FACTOR (DIMENSIONAL) FOR REACTION 1, ETC. 

RHC, MCO = MOLECULAR WEIGHTS FOR CNHM, CO 

OIS = HEAT Of CCMbUSTICN FCR REACTION 1 (CAL/G OF CNHM) 

02S s HEAT OF CCMBUSTICN FOP PEACTICN 2 (CAL/G OF CC ) 

OS = HEAT OF COMBUSTION FOR REACTION 3 ( CAL/G Of CNHM) 

XN = THE NUMBER CF CARBONS IN CNHM 
XM = THE NUMBER CF HYDROGEN IN CNHM 


OS s OIS 4 02S»XN*WCO/WHC 

***** CATALYTIC BED PARAMETERS ***** 

XLS s BED LENGTH (M) 

DS = CHANNEL HYDRAULIC DIAmETER(M) 

SOA = CIRCUMFERENNT lAL LENTH DIVIDED BY CHANNEL AREA,=4/CS (1/M) 
ADAS = RATIO CF OPEN TO CLOSED AREAS 
CSS s heat capacity CF SUBSTRATE (CAL/KG/K) 

R05= DE5ITY OF SUBSTRATE (KG/ M3) 


•READ(5,10) XPHI 
READ (5,10) Tl 
READ (6,10) XFl 
READ (5,10) UlUS 
READ (5,10) XLS 
READ (5,10) DS 
■READ (5,10) ADAS 
PEAD(b,A) UXX, XAS 
READ (5,10) 6 t 
READ (5,10) TMAX 
READ(5,6) N 

RFAD(5,8) NAS 
PFAD(5,8) IT 

RFAD(5,*) C0.C1,C2,C3,C4 
READ(5,8) IPUTIN 
SOA s 4.0/DS 
CN s XnMWCO/WHC 
NPNAS=N4NAS 

****^ OPERATING PARAMETERS 

TI = INITIAL, UPSTREAM (l.U,) TEMPERATURE 

UI = l.U. VELOCITY (M/S) 

XPI s I. U. PRESSURE (ATM.) 

ROIs l.U. gas density (KG/M3) 

YFI S l.U. MASS FRACTION OF HC 

YCI = l.U. MASS FRACTION OF CC 

YOI » l.U. MASS FRACTION OF 02 

XPHI « UPSTREAM HC/02 EQUIVALENCE RATIO 


(K) 


UI » UIUS»(1 ,04AOAS)/AOAS 
XMU r 0.00000145«4SORT(TI)/(1. 04110. 4/TI) 
YFIsWHC*XPH 1/29./((XN4XM/4. )44.764XPH1) 
XLEWl « 1.54 SORT (29. /WHO 
XLEW2=1.5» S0RT(29./WC0) 
XLE123=XLEW1»4.667 
XLE223 sXLEm2»*.667 

XLOU e XLS/Ul 

***** GAS PROPERTY ***** 


CP s SPECIFIC HEAT OF GAS (CAL/KG/K) 

XKI = HEAT CONDUCTIVITY OF GAS AT l.U. (CAL/ M/5EC/K) 
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14200 
1 4300 
14400 
14500 
14b00 
14700 
14750 
14100 
14900 
15000 
15100 
15200 
15300 
15350 
15375 
153*7 
15400 
15500 
15600 
15700 
15100 
15900 
16000 
16100 
16200 
16300 
16400 
16500 
16600 
16700 
16750 
16*00 
16900 
1 7000 
17100 
17150 
17200 
17300 
17400 
17500 
17900 
17950 
18000 
1*100 
1*200 
18250 
18275 
1 8287 
1 6300 
1*400 
1 8500 
1*600 
18700 
1*800 
18900 
19000 
19100 
19200 
19300 
19400 
19500 
1 9600 
19650 
19675 
19700 
19*00 
19900 
19950 
20000 
20100 
20200 
20300 
20400 
20500 
21100 
21200 
21300 
21400 
21500 
21600 
21650 


12-APR-1981 23;31s21.21 PAGE 3 

ALPHA = I.U. GAS THEPHAL UlFEllSIYJTV ( 1'2/SEC) 

GCtGAS CONSTAM EOH GAS H1X7U6F ( JCJJLE/KG/K ) 

XLFWl r LEhJS NUPHEH FCH CNHM IN MIXIUHE 
XLt«2 = LENIE NUPbtH FOP CO IN MIXTURE 


XM = .n03+.00001 3»Tl 
HQl = XHl/GC/TI 
ALPHA = XKl/RCI/CP 
REY = UIUS»DS»RCI/XMU 


TO DESCRIBE chemical REACTION RATE EXPRESSIONS 
REACTION 0 (GAS PHASE ISOTHERMAL PYROLYSIS) 

D(CNHM)/DT S$S (CNHM)/EXP(EO/RI ) 

REACTION 1 { GAS PHASE ) 

D(CNHM)/OT SSS XP * * . 3»1 *5081 ( YHC) ♦ YC2/EXP C E 1 /PT ) 

REACTION 2 (GAS PHASE ) 

D(C0)/D'i $ P4»BTAlMC0»4eTA2»C2**BTA3AT**9TA4*H20»*PTA5/EXPCE2/RT) 

REACTION 3 (SOLIO SURFACE) 

D(CNHM)/DT $$S YHCS/EXP(E3/RT) 

REACTION 4 C50LIC SURFACE) 

D(CC)/DI ($S YCC5/EXP(E4/RI) 

TO EVALUATE KINETIC RATE CONSTANTS. 

BO=CO*XLOU»(HHC/ROI/YFI )♦♦( 1 ,-GAMA ) 

BISS = Cl»S0RT(RHC)/t>C2/SGRT ( 1000. ) 

BIS = B1SS»XP14*.34R0I*S0RT(RCI)*TI»SQRI(YF1 )»YOI 
B1 = XL0U»B1S/RCI/YH 
B235=BTA2+bTA3+BT A5 
C2M=C2A1000.4*(1.-B235) 

B2SSsC2M**iCC»»( 1 .-PTA2)/(*l02**BTA3»»«H204»BTA5) 
b2S=B2SS*XPIM»BlAl*R0IAAB235*TI*4bTA44YFI*»BTA24Y0I*4flTA3 
B2 = B2S4XL0U/RC1/YFI 

B3 = C3»DS/ALPHA/»HC/3 .66 
B4=C4*DS/ALPHA/KC0/3.6b 

01 = QlS/CP/Tl 

02 = 02S/CP/TI 
0 = OS/CP/TI 
TAB=YFl»u+l .0 
TABDMrTABYTI 

E0= ESO/TI/1.987 
El= ESl/Il/1.967 
E2s ES2/TI/1.987 
E3= ES3/II/1.987 
E4= ES4/TI/1.987 
AD2UL = ( ALPHA/LS/DS) /(Ul/XLS) 


MRITE (6 
WRITE(6, 
NRITE (6 
NRITECb, 
WRITE (6 
WRITE (6 
WRITE (6 
WRITE(6, 
WRITECb. 
WRITE (6 
WRITE(6, 
WRITECb, 
WRITE (6 
WRITE (b 
WRITE (6 
WRITE (6 
WRITE (6 
WRITE (6 
WRITE(6, 


,21) N,NAS 
20) TMAX 

,15) YCI, YFl, XPHl 
71) UlUS 
,72) U1 

,17) 01, 02, DXX, DT 
,402)TI ,XPR,REY, AC2UL, ADAS 

409) TAB 

410) TABCM 

,403)ES0,ESi,ES2,ES3,ES4 

47)E0,E1 ,E2,E3,E4 

35)C0,C1 ,C2,C3,C4 

,411) BTA1,BTA2,BTA3,BTA4,BTA5 

, 48) GC, XLEWl, XLEW2, XLS, DS 

,49) ROI, XKl, CP, ALPHA 

,401) XPl 

,43) XLCU 

,44) BIS, B2S 

45)B0,B1 ,B2,B3,B4 




«««•«•»* 

YF(l)sl. 

YE(1)=0. 


UPSTREAM CONDITIONS 
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21700 

21800 

21900 

22000 

22100 

22200 

22250 

22300 

22400 

23000 

23100 

23200 

23300 

23350 

23375 

23387 

23393 

23396 

23397 

23398 
23449 
23474 
23417 
23493 
23500 
23550 
23600 
23700 
23800 
23900 
24000 
24100 
24200 
24300 
24350 
24362 

24374 

24375 
24400 
24500 
24600 
24700 
24600 
24900 
25000 
25100 
25200 
25300 
25325 
25337 

25349 

25350 
25450 
25500 
25550 
25600 
25700 
25100 
26250 
26300 
26350 
26400 
26500 
26600 
26700 
26800 
26900 
27000 
27100 
27200 
27250 
27300 
27400 
27500 
27600 
27700 
27800 
27900 
27950 
28000 
21100 


YCd) = 0.0 
YC(1)=1. 

T(1 ) = 1 . 

R0(1)=1. 

xpcn =1.0 

U(l)=l. 

YTHCCl) = YF(1)+VE(1) 


^*t********:^*f^-*n^************************* 

SIAHT OF INTEGRATION LOOP 
KOUM = 1 

XNUO=.8854SQPT(ACA5)*SORT(XPK*REY) 

ROOT1=0.0 
R0012=l .0 

555 ROOTS. 5»(R00T1+PC0T2) 

IF(RCOT.LE.O.l) YNU=.04 44/R0CI-f 3.46«1 . 34 ♦ROOT 
1F(RC0T.G1 .0.1 ) YNUsO.Oll /ROOT •♦3.66 
YROOT=YNU-XNUO 

IFCAESCYPOOTl.Ll.O.Ol) GO TO 200 
IF(YkOOI.GT.O.O) ROCTl=ROOT 
IF (YkOOT.LE.O.Ol ROOT2 = ROUl 
GO TO 555 
200 CONTINUE 

DHASS = RO( 1 )4U( 1 ) 

INTEGRATION OF GUASI-STEACY GAS PnASE EQUATIONS FOR 
YF (FUEL NASS FRACTION) AND T ( TENPERATURE ) IN DOWNSTREAM 
DIRECTION USING FOURTH-ORDER RUNGE-KUTTA SCHEME 
I IS THE INDEX FCR DOWNSTREAM DISTANCE X 

DO 201 I=1,NPNAS 

IF(l.LE.N) DX=CXX 
IF(I.GT.N) DX=XAS*UXX 

IFCI .GT.N )DMASS = EOC 1 )*U(1 ) *ACAS/ ( AOAS'Rl . ) 

^^^t***’f*********************$*** 

EVALUATE XJH AND XJC AS A FUNCTION CF X 

NONCIMENSIONAL HEAT CCNDUCIIVITY 
XK = (.003't 0.00001 3» TI *1(1)) /XKl 

t$$*$$$$**$f$*$****$$$* 

IF (I.GT.I 1 GO TC 180 
IF(I.EO.l) X=0.0 
GO TC 181 

180 CONTINUE 

IX = I-U 
XI=1X-1 

X = XI*DX 

181 CONTINUE 

XST= X*2.0*AD2UL4XK-fR0CT 

IF (XST.LE.0.1 ) XNU = 0.0444/XST'*3.46 -1 .34*XST 

IFCXST.GT.O.n XNUsO.Oll/XST •♦3.66 

RNKll)=XNU»XK/3,66 

XOH(I) s 4.04XNU ♦AD2UL4XK 

IF(I.GT.N) XJH(I) =0.0 

XJCl(I) « XJh(I) 4XLE123 

XJD2(I) = XJH(I) 4XLE223 

1.32 CONTINUE 
H = 1 • 

DO 120 L=l,4 
IF (N.NE.l) GO TC 102 
YFF= YF (I) 

YEE =YE(I) 

YCC =YC(1) 

TT = T(I ) 

M = M-fl 
CO TO 119 
102 CONTINUE 

IF (H.NE.2) GO TO 103 
YFF=YFCI)-eDX*XKY(1) /2. 

YEE = YE(I)^fDX*XKE(l)/2. 

YCC sYCCI) ■♦DX»XKC(l)/2. 

TT =TCI)^EDX4XKT(1) /2. 
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21200 

28300 

2M00 

28500 

28600 

28650 

21700 

28100 

28900 

29000 

29100 

29200 

29250 

29300 

29400 

29500 

29600 

29700 

29800 

29850 

29175 

29900 

30000 

30025 

30050 

30075 

30100 

30200 

30250 

30300 

30400 

30500 

30600 

30700 

30750 

30775 

30800 

30900 

31000 

31100 

31125 

31150 

31200 

31300 

31400 

31500 

31550 

31600 

31700 

31800 

31900 

32000 

32100 

32200 

32300 

32400 

32500 

32600 

32700 

32800 

32900 

33000 

33100 

33200 

33300 

33400 

33500 

33600 

33700 

33800 

33900 

34000 

34100 

34200 

34300 

34400 

34500 

34600 

34700 

34800 

34900 


M = M+1 
GO TC 119 

103 CONTINUE 
IFCM.NE.3) go 7C 104 

YFF = Vf (I) +DH*XKY(2) /2. 

YEt= YF.Cn'tDX*XFE(2)/2 
YCC =YC(1) ♦DX»XKC(2)/2. 

TT = T(I) + DX’*XKTt2) /2. 

M = M+ 1 
GO TC 119 

104 CONTINUE 

YFF = YF(I) + 0X»XKY(3) 

YEFs YE(I) +DX»XKE(3) 

YCC =YC(I) ♦DX»XKC(3) 

TT=T(I) +DX»XKT(3) 

119 CONTINUE 

YH20 = NASS FRACTION GF H2C 

YH20(I) = .5»XN*8H20*YFI* ( YTHC ( 1 )-YTHC ( I) ) /wHC 

IF (YFf.LT.0.0) YFF=0.0 

WO = BO*(ROm*YFF)**GANA/EXF(EO/TT) 

IF(YEE.LT.O.O) YEL=0.0 

hi = Bl* XP(l)»*.34R0(I)*S0RI(B0m)4T(I)*S(jRT(YEE) ♦YD(I) 

1 /EXP(E1/TT) 

IF(YCC.LT.O.O) YCCsO.O 
SSY0=Y0(I)»*BTA3 
R0235=P0{ 1 )*»B235 

h2=B2»R0235*S5YC*YCC*4BTA2*YH20(I)**BTA5*XPtI)*»BTAl/EXP(t2/IT) 

XKY(L) = -XJCim »(YFF- YF5 ( I ) ) /DN ASS -hO/DNASS 

XKECL)=-XJDim»lYEE-YE5(I))/DNASS+(hO-<il)/DNA5S 

XKC(L) = -XJD2CI) •(YCC-YCStl))/ DKASS ♦CN* w 1 /ON ASS -hJ/DWASS 

XKTCU =-XJH(I)4(TT-T5(I))/ChASS +01AYFI4W1/DNASS+024YFI4«2/DNASS 

120 CONTINUE 

YF(I+1) =YF(1) +CX»(XKY (l)/6. +XKY (2)/3. +XKY (3)/3.+XKY (4)/6.) 
IF (YFU + 1 ) .LI .0.0) YFU + l) sO.O 

YE(I+l)=YE(l)+DX4CXKE(l)/6.+XFE(2)/3.+XFE(3)/3.+XKE(4)/6.) 

IFCYEd + D.LT.O.O) YE(1 + 1) = 0,0 

YC(I+1) =YC(1) 4CXACXKC (l)/6, +XKC (2)/3. +XKC (3)/3.+XKC (4)/6.) 
IF (YC(l+n .LT.0.0) YC(I+1) =0.0 

T (I + n =T (1) +LX»(XKT m/6. +XKT (2)/3. +XKT (3)/3.+XM (4)/tj.) 
IF (T ( I + l ) .LT . .0 ) GO TC 601 
YTHC Cl+n = YF(I + l) + YE( I + n 
ADD=YC( 1+1 )*YFI»»02/YCI/hCC/2. 

Y0(I + 1) =1.0-XPhl + (1.0-YfCI + n) + ADC 
XP(1+1)=XP(1) 

ROtl + 1 )=XP(I + n /T(I + 1 ) 

U(l+1) =DNASS/Tt]+1) 

EFF( n = l .-YThCCI )-YCCl)»G2/0 
201 CONTINUE 

END OF ONE RUNGE-KUTTA INTEGRATION CYCLE FOR OU AS I -STF. A D Y 
GAS PHASE EQUATIONS 


««*««==*=**=«***+»«A**«**«**«*«4*«+**4*+A*» 

INTEGRATION OF UNSTEADY SOLID NEAT TRANSFER EQUATION 
FOR IS C SOLID TEMPERATURE ) IN TIME USING FOURTH-ORDER 
RUNGE-KUTTA SCHEME 


DO 301 J *1,N 
MJ=1 

DO 320 LJ=1,4 

IF (MJ.NE.l) GO ID 302 
TSS =TS(J) 

MJsMJ+1 
GO TO 319 

302 CONTINUE 
IF(MJ.NE.2J GO TO 303 

TSS = TS(J) + DT*XKTSm/2, 

MJsMJ+1 

GO TO 319 

303 CONTINUE 
1F(MJ,NE.3) go to 304 
TSS=IS(J)+DT»XKTS(2) / 2 . 
MJ=MJ+1 

GO TO 319 

304 CONTINUE 

TSS = TS(JJ + DI4 XKTE(3) 
319 CONTINUE 
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34950 

35000 

35100 

35200 

35300 

35400 

35450 

35500 

35600 

35700 

35800 

35900 

36000 

36100 

36200 

36300 

36400 

36500 

36600 

36700 

36SOO 

36900 

36950 

37000 

37100 

37200 

37300 

37400 

37500 

37600 

37650 

37675 

37687 

37699 

37700 
37800 
37900 
31000 
38100 
3*200 
38300 
38400 
38500 
31543 
38586 
3*593 
38594 
3*596 
3*599 
3*624 
3*650 
3*700 
3*800 
3*900 
39000 
39100 
39200 
39300 
39400 
39500 
39600 


YfSES=VFS(J)4YES(J) 

+B3AYFI»Q»PO(J)AYFSES/EXP(E3/7SS) 

1 4B4*YFI»Q2*R0(J)»YC5(J)/FXP(E4/TSS) 

320 CONIINUE 

+DT*(XI<TS( 1 )/b. + XI<TS(2)/J. + XKlS(3)/3.4XKTS(4)/b. ) 
/(I .483ARC(J)/XtE123/RNKlJ)/EXP(E3/TS(J) n 
YES(J)=YE(J)/(1.4B34RC(J)/XLE123/RNK(J)/EXF(E3/15CJ))) 

301 CONTINUE /(1.+B4»R0(J)/XLE223/RNK(J)/EXP(E4/TS(J))) 

LOOP FOR MARCHING FORWARD IN TIME 

KOUM bKOUNT+I 

KOU sKOUNT-1 

XKOUNT =KOUNT 

TIME=(XK0UNT-1.)»DT 

IF (TIME.GT. TMAX) GO 10 999 

TIMER=TMAX-DT 

IFCTIME.lt. TIMER) GO TO 200 
WRI1EC6,22) TIME ,KOUNT 
WHITE (6,13) 

WR1IE(6,24) (T CI),I=1,NPNAS,IT) 

WRITE (6,9) 

WRITE(6,24) (TS (I), 1 = 1, N, IT) 
wRITE(6,14) 

WRITE(6,24) (YF(I),I=1,NPNAS,IT) 

WRIT£(6,412) 

WRITE(b,24) (YE(I),1=1 ,NPNAS,II) 

WRIT£(6,413) 

WRITE (6, 24) (YTHC(I),I=1,NPNAS,IT) 

WPITt(6,41) 

WRITE(6,24) 

WRITE(6,11) 

WRITE(6,24) 

WRITEC6,42) 

WRITE(6,24) 

WR1TE(6,26) 

WRITE(6,24) (XJH(I),I=1,N,IT) 

WRITE(6,73) 

WRITE(6,24) 

WR1TE(6,23) 

901 CONTINUE 

IF (IPUTIN.EG. 1 ) GO TO 200 
0PEN(UMT=1 1 ,TYPE='NEW',NAME= 'SOLID') 

WRITE (!!,♦) A 
CLOSE (UNITsll) 

GO TO 200 

END OF WHOLE INTEGRATION LOOP 


(YC(I),I=1,NPNAS,IT) 

(YFE(I),I=1,N,IT) 

(YCECI),I=1,N,IT) 


(EFf (1 ) ,1=1 ,NFNAS, IT) 


601 CONTINUE 

WRITE (6,404) 
WRITE (6,405) 
999 CONTINUE 
STOP 
END 
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100 
200 
300 
400 
600 
700 
• 00 
900 
1100 
1150 
1175 
1200 
1400 
1500 
1600 
1700 
1(00 
1900 
1950 
2000 
2100 


BLOCK DATA 

dimension 7S( 300 ), YFS( 300), vest 30 U) , yes ( 300) 

COMf*CN /AClI/tSO,tSl,ES2,tS3,ES4 
COMMON /MULFO/ mHC, WCC, KU2, KH2U 
COMMCN /IHl/XN, XH, OtS, 025 
COMMON /TH2/ CF, GC, XPR 
COMMON /UPST/ YCl, KOI 
COMMON /INI/TS,YFS,KCS,VES 
COMMON /B12345/B1A1 ,BlA2,BTA3,BTA4,eiA5 
COMMON /GAMMA/GAMA 

DATA ESO.ESl ,ES2,ES3,ES4/34000., 24000, ,40000. , 10000. ,17800./ 
DATA MhC,MC0,i«C2,MH20/l3«.,28., 32.,ll./ 

DATA X N,XM, 015, 025/ 10.,1«. ,5357 000. ,2415000./ 

DATA CF,GC,XPP/2e5. ,2.*312E-3,0.7/ 

DATA VCI ,KOI/0.,0.23/ 

DATA E(TAl,BTA2,E:TA3,eTA4,BTA5/0.,l., 0.25,0., 0.5/ 

DATA GAMA/0,5/ 

DATA T5,KF5,KCS, YES/ 30041., 300*0., 300*0., 300*0./ 

END 
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Example calculation 


r 


S RUM CRLY1.exe; 61 
.165 
1 0 0 0 . 

15. 

. 1 

. 00138£ 

13. £9 
.01 1 . 

. 1 

8 . 001 
51 
1 0 0 

c 

1.6E6 1.8E5 .54E14 1 . 5E4 1 . E5 
1 

ND. DF GRID PDIMT = 51 







100 





4 

71 

xpi 

UlUS 

)(LS 

PS 

AO AS 

3>XX ^ XAS 

DT 

TMAX 

NAS 

IT 

C0,CI,C2,C3.C4 

IP'CTTtN 


TMFIX = 8. 00100040 

'.■□aiPST» INI > = 0. £30000 YF OJPST.INI:) = 0.0113489 EOUI RflTID = 0.165 

00000 ? 


BPPRDftCH VEL = 15.0000 

CHRNNEL VEL = 16. 1£867 


01= 18.7964916 Q£= 8.4736843 DX= 0.0100000 DT= 0. 

1 0 0 0 0 0 0 


Til. PRj REY, PD£UL, ftDfiS = 1000. 00000 0.70000 589. 01978 0.17199 13 

. £9000 


FiDIPiBflTIC FLftME TEMP = 1.4 0844 0 

FIDIflB FLRME TEMP <tO = 1408.44019 
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ES:-S = 
OCiE + 05 


0. 34 0 0Cnji:M;iE+0? 

0. 17800i!iOCiE-*-0“ 


0 . £4 0 Cl 0 0 0 OE + 0? 


0. 4CiCiCiiKiCiCiE + 0? 


Cl . 1 Cl Cl Cl Cl Cl 




17.111E23 1£. 078T.1 Cl £Ci. 1 SCiSSCi 


Ci3c71c' 95fii£’c'& 


CjS = Cl. 16 Cl 000 Cl Cl 00E+ 07 0. ISCiCiCiOCiOCiCiE + OE 0. 5399999866E+14 0.1500 

0 0 0 0 0 0 E + 0 5 0. 1 0 0 0 0 0 0 0 0 0 E + 06 


BiETEl 1-5 = 0.0000 1.0000 0.2500 0. 0000 0.5000 


GCn XLEWli. XLEW£» XLSu DS = 0. 00E831 0.687624 1.526551 0.1 00000 0 

. 001 382 

801. XKI. CP. flLPHP = 1.059621 0.016000 285. 000000 CI.52982E-04 


XP I = 3. 000 0 0 

L- U = 0. 0062001 

BUS. E2S. = 0. 7764775 OE-t- 05 0. 261 7391 0E-*-09 * 

pO-4 = 0. 1 0626966E407 0. 40033785E+05 0. 13494798E-*-09 0. 77466638E + 03 0. 25453324E+0 


TIt1E= 8. 00000000 KDUMT= 81 


TEMPERATURE CT> 

0. 1 OOOOE+Ol 0.1CI933E + 01 0. 1 1328E+01 0. 1 1648E401 0. 11927E+01 0. 12180E + 01 0. 12416E 
+ 01 0. 12640E+CI1 0.12855E+01 0. 13063E+01 

0.13262E+01 0.13387E+01 0.13506E+01 0. 13629E+01 0.13739E+01 0.13806E+01 0.13853E 
+01 0. 13891 E+ 01 0. 13923E+01 0.13948E+01 

Ci.r3969E + 01 0.13985E + 01 0.13998E+01 0. 14008E+01 0. 14016E+01 0. 14022E+01 0. 14027E 

+01 0.14031E+01 0.14034E+01 0. 14036E+01 

0.14038E+01 

SOLID TEMPERATURE <.TS> 

0. 1 1286E+01 0.13193E+01 0.13325E+01 0.13421E+01 0.13502E+01 0.13576E+01 0.13646E 
+ 01 0.13713E+01 0. r5779E + 01 0. 1384£E + 01 
0. 13901E+01 

□RI6IMAL CMHM CYP> 

O.lOOOOE+01 0.76186E+00 0.65814E+00 0.57097E+00 0.49212E+00 0.41866E+00 0.34867E 
+00 0.28172E+00 0.21825E+00 0. 15913E+00 

0.10576E+00 0.65733E-01 0.33518E-01 0.10262E-01 0, 16875E-03 0. OOOOOE+00 0. OOOOOE 
+ 00 0. OOOOOE + 00 0. OOOCiOE+OO 0. OOOOOE+OO 

0. OOCiOOE+OO 0. OOOOOE+OO 0. OOOOOE+OO 0. OOOOOE+OO 0. OOOOOE+OO 0. OOOOOE+OO 0. OOOOOE 
+00 0. OOOOOE+OO 0. OOOOOE+OO 0. OOOOOE+OO 
0. OOOOOE+OO 
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PYRDLY2EH CNHM CYE' 

0. OOOOOE + 00 0. 37Ci££E-0£' 0.91331E-0E' 0. 15t.43E-Cil 0. c’E'801E-Cil 0.3Cilc7E-01 ij. 3703t.E 
-01 0.4£908E-01 0.47043E-01 0. 48€.£7E-01 

0.48778E-01 0.4431£E-01 0.36474E-01 0. ££91 OE-01 0.59888E-0£ 0. OOOOOE+00 0. OOOOOE 
+ 0 0 0 . 0 0 0 0 0E+ 0 0 0 . 0 0 0 0 0E+ 0 0 0 . 0 0 0 0 0E+ 0 0 

0. OOOOOE+OO 0. OOOOOE+00 0. OOOOOE+00 0. OOOOOE+00 0. OOOOOE+00 0. OOOOOE+00 0. OOOOOE 
+00 0. OOOOOE+00 0. OOOOOE+00 0. OOOOOE+00 
0. OOOOOE +00 

TDTPlL CNHM (YTHO 

O.lOOOOE+01 0. 785'56E+00 0.867£7E+00 0.5886£E+00 0.51493E+00 0.44878E+00 0.38571E 
+00 0.3£483E+00 0.£85£9E+00 0.£0778E+00 

0.15£54E+00 0.11004E+00 0.6999£E-01 0.3317£E-01 0.81574E-0£ 0. OOOOOE+00 0. OOOOOE 
+00 O.OOOOOE+00 O.OOOOOE+00 0. OOOOOE+00 

0. OOOOOE+00 0. OOOOOE+00 0. OOOOOE+00 0. OOOOOE+00 0. OOOOOE+00 0. OOOOOE+00 0. OOOOOE 
+00 0. OOOOOE+00 0. OOOOOE+00 0. OOOOOE+00 
0. OOOOOE+00 

Yca = 

0. OOOOOE + 00 0. 16567E-0E' 0.853£4E-0£ 0. 15098E-01 0.£7740E-01 0.44585E-01 0.6?'i£77E 
-01 0.89093E-01 0.114S7E+00 0.14085E+00 

0. 16454E+00 0.£1374E + 00 0.£6086E+00 0.£8957E+00 0.£8939E+00 0.£4845E,+ 00 0. 198.388 
+00 0.15831E+00 0. l£54£E+00 0.98789E-01 

0.77393E-01 0.80398E-01 0.48983E-01 0.38453E-01 0.£8££8E-01 0.£18£1E-01 0.18848E 
-01 0.1£998E-01 0.10017E-01 0.77188E-0£ 

0.59417E-0£ 

SLIRFHCE HC M8SS FRRCTIDN (YFS> 

0.59588E+00 0.555£7E-01 0.41£78E-01 0. 33799E-01 0.£8303E-01 0.£3840E-01 0.£007?;E 
-01 0.18398E-01 0.1£835E-01 0.94545E-0£ 

0.63485E-0£ 

YCD':SURFftCE> = 

O.OOOOOE+00 0.1££48E-03 0.404£4E-03 0.88451E-03 0.15178E-0£ 0.£3790E-0£ 0.34733E 
— 0£ 0-47£8£E— 0£ 0-80618E— Oc 0»744£OE— 0£ 

0.88843E-0£ 


XJH 

0.4£719E+0£ 0.3989£E+Ul 0.33666E+01 0.31814E+01 0.308£3E+01 0,30479E+01 0.308£7E 
+01 0.31189E+01 0.31558E+01 0.319££E+01 
0. 3:££83E+01 

EFF I C I ENC't <:CftRIiDN-BFtLFlNCE> 

0. OOOOOE+00 0.83405E+00 0.33119E+00 0.40983E+00 0.47854E+00 0.54072E+00 0.59893E 
+00 0.65439E+00 0.70766E+00 0.75908E+00 

0.8087£E+00 0.83983E+00 0.88884E+00 0.89885E+00 0.9£571E+00 0.94197E+00 0.95330E 
+00 0.98£73E+00 0.97047E+00 0.97674E+00 

0.98178E+00 0.98578E+00 0.98894E+00 0.99143E+00 0.99335E+00 0.99488E+00 0.99603E 

+00 0.99694E+00 0.99764E+00 0.99818E+00 

0.99860E+00 


83 



1 Report No 

NASA CR-16532/. 

2 Government Accaision No 

4 Title and Subtitle 


Transient Catalytic Combustor Model 

7. Author{s) 


James S. T'ien 


9. Performing Organization Name end Address 


Case Western Reserve University 
2040 Adalbert Road 
Cleveland, Ohio 44106 


12. Sponsoring Agency Name and Address 


U.S. Department of Energy 
Office of Transportation Programs 
Washington, D. C. 20545 


3 Recipient's Catalog No 


t Retxjri Date 
Mav 1981 


6 Pefforrnmg Ot 9 ani?ation Cooc 


8 Performing Organization RepoM No 


10. Work Uni! No 


11. Contract or Grant No. 
NSG - 3230 


13. Type of Report and Period Covered 
Contractor Report 


14. Sponsoring Agency Code 
DOE/NASA/3230-l/l 


15 Supplementary Notes 

Final Report, prepared under Interagency Agreement DEAI01-77CS51040. Project Manager, D. Anderson, 
Aerothermodynamics and Fuels Division, NASA Lewis Research Center, Cleveland, Ohio 44135.' 


16 Abstract 

A quasisteady-gas-phase and thermally-thin-substrate model is used to analyze the transient 
behavior of catalytic monolith combustors in fuel-lean operation. The combustor response delay 
is due to the substrate thermal inertia. Fast response is found to be favored by thin substrate, 
short catalytic bed length, high combustor inlet and final temperatures and, in most cases, 
small gas channel diameters. The calculated gas and substrate temperature time history at 
different axial positions provides an understanding of how the catalytic combustor responds 
to an upstream condition change. The computed results also suggest that the gas residence times 
in the catalytic bed and in the after-bed space are correlatable with the nondimensional combustor 
response time. The model can also perform steady-state combustion calculations; and the computed 
steady-state emission characteristics show good agreement with available experimental data in 
the range of parameters covered. From combined steady and unsteady considerations, it is possible 
to choose a catalytic combustor design for automotive gas turbine engine which has reasonably 
fast response ( < 1 second) and can satisfy the emission goals in an acceptable total combustor 
length. 


17. Key Words (Suigested by Author|s)l 

Catalytic combustion 
Transient combustion 
Combustor modeling 


IS. Distribution Statement 

Unclassified - unlimited 
Star Category 44 
DOE Category UC-96 


19. Security OiKif. (of this report) 
Unclassified 


20. Security Qessif. (of this page) 

21. No. of Pages 

22. ^icc* 

Unclassified 




* For sale by the National Technical Information Service, SprinElield. Virginia 22161 


NASA-C-16I (Rev. 10-75) 


84 























CENTER 


II I ill ill III I 

lllilililll 


3 1176 00507 3466 


» 


